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The present thesis is an attempt, to evaluate the 
performance of three different semi-empirical MO formalisms 
vis-a-vis the rigorous evaluation of electric field gradients 
(BPG) by ab initio quantum mechanical -procedures and experiments. 
The formalisms selected for study are the CNDO/2 and INDO methods 
of Pople et al. , and a method due to Nanda and Narasimhan at the 
INDO level which uses a basis of Orthogonalised Atomic Orbitals 
referred to as the NN-INDO method. This method achieves rota- 
tional invariance while differentiating the s and p orbitals, 
and the p orbitals among themselves. EPG’s at the sites of the 
quadrupolar nuclei D, Li, B, N and 0 in some selected molecules 
are studied. The EPG's are evaluated rigorously using calculated 
MO wave functions in these formalisms for the molecules under 
study, and the results are analysed in the light of other existing 
calculations and experiments to obtain an insight into the physico- 
chemical phenomena involved. 



The thesis consists of five chapters and two appendices. 
References to the literature a,ppear at the end of each chapter. 

Chapter I forms an introduction to the work desfrihed in 
the thesis. The nature of the quadrupolar interaction is out- 
lined and ab initio as well as empirical and semi-empirical 
methods of evaluating EPG's are critically evaluated. 

Chapter II outlines the semi-empirical formalisms chosen, 
and describes the gaussian transform method used for calculating : 
the integrals required for evaluating the EPG-'s. 

Chapter III contains a study of the EFG's at the sites of 
quadrupolar nuclei in the linear molecules lithium hydride, 
hydrogen fluoride, carbon monoxide, nitrogen, hydrogen cyanide, 
cyanogen and acetylene. It is found that the deuterium field 
gradients obtained by the RN-INDO method are consistently higher ; 
than that obtained by the other two methods, and this is accoun- ^ 
ted for on the basis of the Burns exponents that are used in the i 
M-IRDO method. -j 

Chapter IV deals with the calculation of EPG- parameters I 
in some non-linear molecules. The molecules chosen.: for study I 
are methane, ammonia, water, hydrazine and formaldehyde and i 

the cyclic molecules pyrrole and furan. Both the field gradient ' 
and the asymmetry parameter are evaluated at the sites of the j 
quadrupolar nuclei# j 

Chapter V contains our results using the three different 
semi-empirical formalisms for the various quadrupolar nuclear 



sites in the molecules ethylene, methyl cyanide, methylamine, 
formamide, hydrogen peroxide, nitrous acid, diborane, borazine, 
boron trifluoride, nitrogen difluride, fluorine monoxide and 
fluorine dioxide. Some interesting qualitative correlations 
involving chemical concepts like hybridisation, electronegativity, i 
etc. are observed and reported along with the results for the 
field gradient parameter. 

The thesis concludes with a few critical remarks on the use 
of semi-empirical methods for evaluation of EPGr's, 

Appendix I of the thesis gives a listing of a computer program; 
used in this work for the evalua.tion of the diagonal component of 
the EPfi, Appendix II lists the conversion factors used for the 
conversion of experimental nuclear quadrupole coupling data into 
EPG’s. 
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I IVTRO RUCTION 
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1.1 The Mature of the Quadrupolar Interaction 

A wide class of interesting phenomena have their origins 
in the interaction of the nuclear charge distribution with the 
surrounding electron distribution. Among these may be listed 
the hyperfine structure and isotope shifts in high resolution 
spectroscopy, quadrupole hyperfine structure in rotational 
spectra of molecules, quadrupole splitting in Mossbauer lines., 
and the phenomenon of Nuclear Quadrupole Resonance (NQR). 

Consider a nuclear charge distribution given as p(l'). The 
electron distribution surrounding the nucleus generates an ele- 
ctric potential V(r). The interaction energy is given by 

E = Jp(r)’V(r)dr (1.1) 

Expanding the potential VCr) in a Taylor series about the origin, 
we get 

E = V(o) I p(r) dr + I (|~ ) x^^^ P(r) dr 



higher order terms (1.2) 

The first term gives the Coulomb 'monopole' energy which is 
independent of nuclear orientation. The second term, namely 
the electric dipole moment contribution, vanishes due to nuclear 
parity^. In fact, it can be shown that, all nuclear electric 
moments of odd parity vanish. Hence, to fourth order, the inte- 
raction energy consists of the monopole and the quadrupolar terms 
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This quadrupole terin may be written as 


hj - 3 «lo) '’(r) + 

1,3 

F J. hj ^ '5*' (1-5) 

^ > 0 

The last term is rotationally invariant and hence is not of 
interest to us here. Defining 

Qj_j = I (3Xj_ - r^ 5j_^ ) p(r) dr (1.4) 


We got 


E 


Q 


6 hd 




'E 


( 1 . 5 ) 


This expresses the quadrupolar interaction energy as the pro- 
duct of two dyadics, the nuclear quadrupole moment 

Q = I ( 3 r r - r^ 1 ) p ( r ) dr 


and th^^ gradient of the electric field due to the external 
electron charge distribution having components 


2 

3 V 

(VB),. , = - 7 , . = - (r- — rr ) ,. 
' 13 13 axj^ ax^ -^o 


where the subscript 0 on 


2 
3 V 


9x. 3x. 

4- J 


( 1 . 6 ) 

denotes that this quantity 


is evaluated at the nuclear centre. 


Since the field B is produced by charges external to the 

nucleus, the field gradient tensor v E is traceless, i.e., 

2 

V satisfies the laplaice's equation 7 7 = 0, further, since 
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2 

3 ¥ 


2 

3 ¥ 


3x^ 


3 X . 3 X. 

V 


, it is symmetric. Hence the electric 


field gradient (BPG) tensor has five independent components. 
Also, we can choose the coordinate system such that the off- 
diagonal elements vanish. Such a coordinate system is referred 
to as the principal axis system for the EPG tensor. When the 


field gradient has cubic symmetry, ¥ 

Lanlace's equation ¥ + ¥ + ¥ 

^ XX yy zz 


= ¥ = ¥ . 
XX yy zz 

0; hence, ¥ 


Now, by 


XX yy 

In other words, the quadrupolar interaction vanishes. 


¥ 


zz 


0^ 


The nuclear quadrupcle moment tenser Q is also traceless 
and symmetric, as can be seen from eqn. (1.4). It can be shown 
that the Q tensor has effective cylindrical symmetry and hence 
= Qyy, and =0, if i j . Since I Qj_j_ = 0, we have 
= Qyy = - ^ ^zz' Thus, a single parameter is sufficient to 
specify the nuclear quadrupole moment, and, semiclassic ally we 
can write the interaction energy as 




i 

6 


Q.. (V 




ZZ 


1 


(V 


XX 


¥ ) ) 

yy 


(1.7) 


To obtain the quantum mechanical Hamiltonian from the classical 
energy given by (1.5) one uses the quantum mechanical operator 
for Qj|_j • This gives the Quadrupolar Hamiltonian as 


H 



I 7.. Q, 


(op) 


i , 3 


10 "ij 


( 1 . 8 ) 


To obtain the quadrupolar interaction energy, we need to evaluate 
the matrix elements of H^ between the nuclear states |I,mj^ , 
where I is the spin angular momentum quantum number of the 



5 


nucleus, and is the z-component of the spin which has (21+1) 
•values. Since the energy difference hetween the ground, and 
excited states of nuclei is far greater than Eq, only the ground 
state of the nucleus need be considered. Hence, we need to 
evaluate only such matrix elements as 



2 

Since by the Wigner-Eckart theorem the matrix elements of all 
traceless and symmetric second rank tensors are proportional, 
we may express the matrix element above in terms of the matrix 
element of another operator with the same symmetry as > in 
particular we can write 

I Q°P I 1.4) = C /l.mjl |(I^ I. + I. I.) - 6^^ 1,4) 

(1.9) 

where C is a constant independent of mj, m^, i and j. To 
evaluate C we evaluate the matrix element in eqn. (1.9) for 
mj = = I and i = j = z. 

Defining the scalar quantity Q such that 

eQ = </l,I 1Q°P I 1,1^ (1.10) 

now gives from eqn. (1.9) 

eQ = C <J,I I 31^ - I 1,1^ 


= CI(2I - 1 ) 


( 1 . 11 ) 
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For I values less than 1, eQ vanishes. Hence only nuclei with 
I 1 need be considered here. From eqn.. ( I, 1 1 ) we obtain, 


jQ ,/l(2I - 1 ) 


( 1 . 12 ) 


The quadrupolar Hamiltonian may be now written as 

= TT7o^~7T y V. . (| (I. I. + I. I. ) - 6. . 1“^) 
Q 61(21-1) / ^ ij '2 ' 1 3 3 1 ^ 


1, J 




(1.13) 


Since, the energy is independent of the coordinate system, we 
can now transform to the principal axis - system of the EFG- ten- 
sor. Since, in this system 7 . . =0, if i 3 , and as \ V., = 0, 

XJ XX 

we get 

y - F) + (7^^ - Vyy)(4 - I^)] 


Introducing the qusxntities 

eq = V, 


zz 

V 


and r\ 


- 7 

XX yy 


7 


(1.14) 

(1.15) 

(1.16) 


zz 


where by convention the axes are chosen such that 

1 7 I > 1 7 I > 17 I 

' -zz' ' yy ' I XX I 


gives us the quadrupolar Hamiltonian as 

= {e^q Q/4I(2I-1) } [(34 " + ^(4 - 4^^ (1.17) 


By definition, eq is the component of the EFG tensor having 
the largest magnitude, and ti is the 'asymmetry parameter’ having 
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rtamerical values between 0 and 1. ri represents the departure 

of the EPG tensor from axial symmetiy. These two quantities, 

2 

VIZ., e qQ and p, together with the three Euler angles needed 
to specify the orientation of the principal axes, describe the 
EPG tensor completely. 

2 

The quantity e qQ is known as the ' Quadrupole Coupling 
Constant (QCC), QCC and p are experimentally measurable by 
various techniques. Prom the chemists' point of view these 
quantities are of interest, since they reflect the nature of 
the electron distribution around the nuclei. 

1.2 Experimental Determination of Quadrupole Coupling Constants 
and p ' s 

Several experimental methods are available for the deter- 
mination of Nuclear Quadrupole Coupling Constants (NQCC's) and 
p's. Por studies in gas phase, the techniques of microwave 
spectroscopy, molecular beam magnetic resonance and molecular 
beam electric resonance methods are commonly used. Por studies 
in liquid phase. Nuclear Magnetic Resonance (NMR) relaxation 
data are of some utility. Por studies in solid state, NMR, 
Electron Spin Resonance (BSR), Nuclear Quadrupole Resonance (NQR) 
and Mfossbauer effect techniques are commonly employed, Imongst 

these, the NQR technique has been the most widely used, since the 

3 

main interaction probed is the quadrupolar interaction, lucken 
has given an exhaustive review of the various experimental 
methods of determination of NQCG ' s and p's. 
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The HQCC's and ti's are very sensitive to the state of 
aggregation of the molecules studied. Since most of the mole- 
cular quantum chemical calculations are done on isolated 
molecules, it is preferable to compare the theoretically pre- 
dicted QCC ' s with gas phase measurements. From pure NQR studies 
in solids, it has been clearly established that the quadrupole 
interaction energies can vary considerably depending on inter- 
molecular interactions'^. The measured QCC ' s are fairly accu- 
rate, since the techniques of radio frequency spectroscopy 
are generally employed. In the past decade, tremendous 
strides-^ have been made in the precise determination of QCC ' s 
and Ti's of several nuclei, in a variety of molecules. The 
quantitative interpretation of these data is a matter of great 
interest to the theoretical chemist. In the next section, we 
shall highlight the theoretical aspects of the calculation 
of BPG’s in molecules. 

1.3 Ah Initio Evaluation of Electric Field Gradients in 
Molecules 

The EFG tenser is determined by the electronic and nuclear 
contributions. The classical expressions for the EFG tensor 
components may be obtained treating the nuclei (other than the 
one at which the EFG is to be evaluated and which will be de- 
noted by subscript A) and the electrons as point charges. The 
potential at the nucleus A is given by 

? = I ^ 
p ^AP 


(1.18) 
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where Z is the charge on the particle p and r^^ is the dis- 
tance between the nucleus A and the particle p. The corapo- 
nents of the EFG tensor may be obtained readily from (1. 18)' as 

Xi. 


V. . 
10 


? 


3 


= I Z 

P ® ■ ^Ap 


(1.19) 


Quantum mechanically the quantities V. . are then given by the 

/ \ 10 

expectation value of as given by (1.19). Within the 

10 

Born-Oppenheimer approximation we can treat the nuclei as 
stationary and write 


V. . 

10 


e I Z. 

B;^A 


^3AB - hi 


■AB 


/vj I ( ^^iAS ^0A3 ~ ‘^i.i ^AS 


) Iv) 


( 1 . 20 ) 


Here B runs over all the nuclei other than A and S runs over 
all the electrons in the system. In particular, for we 


have 


eq 


zz 


3 cos e 

e I Z ( - 1 ^“ 


1 


K'i'i n 


■3 cos 0 


AS 


1 


S 




•"AS 


) 


( 1 . 21 ) 


nucl. elec. 

'V in these equations is the total electronic wave function 
for the system. Several choices are available for the molecular 
wave function in eqn. (1.21). For Instance, we may have a 
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valence-bond or a molecular orbital or a correlated etc. type 
of wave function. If the molecular wave function is adequately 
described by a single Slater determinant of molecular orbitals 
(MO), i.e., considering the closed-shell case, we have 

w ^ U a(l) c() |3(2) .... 4 , p(2n)I (1.22) 

where the (J)' s are individual MO's and 2n is the ' total number 
of electrons. Bach of the c{)' s may be expanded in terms of a 
set of atomic orbitals (AO) Xjj- 


{ °ik \ 


(1.25) 


This is the well-known ICAO MO approximation. With this 
approximation eqn. (1.21) reduces to 




nucl. 

f I. Ze { 




ry 2 ^ - OOC 

^ ^AB - 1 .III 

BM B ' rh ^ i k 1 


' AB 


(1.24) 


where 


3 cos^ e - 1 
^Ikl^A “ i" ^ ^ 


'k ^1 


(1.25) 


■A 


Here the index i runs over all occupied MO's and k and 1 run 
over AO's. 


Defining P, the charge-and bond-order matrix, as 
occ 

(B)kl = J, 2 «ik =11 


(1.26) 



and considering the quantities as forming a matrix q, eqn. 

(1.25) may also be written as 


eq 


zz 


A 


rM ® 


3 cos 0 


AB 


1 


-) - € Tr(P q^) (1.27) 


"AB 


Similar expressions are obtained for the other EPG- components 
as well. 


An exact calculation of the EFG- components for a molecular 
system in the ICAO-MO framework requires the evaluation of in- 
tegrals of the type 

IXk q^^ dx =JX^ X^ dx (1.28) 

^A 


Since q^^ is a one-electron operator these integrals may be of | 
one-, two-, or three-centre variety depending on the particular 
AO's involved in each case. The task of evaluating these 
integrals increases in going from one- to two- to three-centre ; 
type. While the evaluation of one-centre type integrals with ■ i 

Slater type orbitals (STO) is fairly straight forward, the 1 

evaluation of two- and three-centre integrals is rather difficult! 

-3 ^ 

Due to the r operator form, closed expressions for two- and j 

three-centre integrals over STO’s are not generally available I 

in the literature. In fact, complete evaluation of all integrals! 

7 i 

was not routinely feasible till the work of Kern and Karplus 
appeared. Prior to this, studies were made mainly on diatomic 
molecules, where only one- and two-centre integrals need be 


evaluated. 



Q 

With reference to ab initio calculations , one of the 
most important factors to be considered is the nature of the 
basis set used. The most common basis set for molecular cal- 
culations are the nodeless STO's. These, however, are not 
altogether satisfactory in simulating the Hartree-Pock AO's 
when used singly. Therefore, extended basis sets are used. 

Prom the point of view of computational simplicity the gaussian 

Q 

wave functions are to be preferred. Many authors have therefore 
expanded STO's in terms of gaussians and proceeded with the 
evaluation of the integrals. The gaussian functions suffer from 
the fact that they have poor cusp behaviour. In order to im- 
prove upon the quality of the Hartree-Pock molecular wave fun- 
ctions, configuration interaction (01) methods have been employed 
It must, however, be pointed out here that the number of inte- 
grals to be evaluated becomes very large with such wave functions 
and computations are therefore generally limited to small mole- 
cules. With particular reference to EPG- calculations using ab 
initio methods we may refer the reader tc an excellent recent 
review by Moccia and Zandomeneghi^ The theoretically calcula- 
ted PPG values in molecules are very sensitive to changes in the 
geometry. Detailed studies have shown that , the 'frozen nucleus' 
calculations have to be supplemented by vibrational averaging 
over the nuclear motion, in order to obtain good agreement with 
experimental results. Such calculations typify the sophistica- 
tion to which the ab initio methods have been developed for the 
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calculation of EFG-'s. With the availability of such precise 
calculations of SPG’s one can even obtain the value of the 
nuclear quadrupole moment from the experimental QCC's. This 
is particularly true of the quadrupole moment of deuterium, 
since very accurate hydrogen molecule wave functions are avai- 
lab le . 

1-4 Bmnirioal and Semi-Empirical Estimates of EFG's in Molecules 

Since, calculations of the EEG- parameters by ab initio 
methods were not feasible in the earlier years, historically', 
empirical and semi-empirical determination of the BEG parameters 
in terms of qualitative valence concepts came first. The pioneer- 
ing work in this direction was that of Townes and Dailey 


Referring to eqn. (1.25) the electronic contribution to the 
BEG at the nucleus of interest, A may be written as 
, ^ occ on A - 


elec . ,A 


-2 I I 

i k 


'ik ^'^kk^A 


2 I T 


on A on 

I ^ ^ °ik ^il^'lkl^A 

B/^A C 

il ^^ll^A + 2 I I I °il'"^im^^lm^A 


i 1 


(1.29) 


m 


i refers to the occupied MO's. The index k refers to orbitals 
on atom A while l,m refer to orbitals centred on other atoms. 
This makes explicit the ' centeredness' of the various integrals 
involved. 



u 


The Townes-Dailey approximation may he now summarised as 
follows 

bM 

1 . The terms ^ '[ assumed to cancel with 

1 1 

the nuclear contributions from atom B. 


2. The quantities (q]^ 2 .)j;, neglected in the simpler versions 

of the theory. This amounts to neglect of contributions of 
overlap densities involving, one orbital on the quadrupolar 
nucleus with the other on another atom. However, these 
contributions were estimated by Townes and Dailey to account 
for only 5 percent of the total BPG. 


3. The three-centre integrals s-re neglected. 

As a consequence of these assumptions, we are now left with 


only the first term on the r.h.s. of eqn. (1.29). In the one- 
centre terras, we have integreils such as I du | 

where k and k' are orbitals centred on A. By symmetry, integrals! 

and 

Xjj.* respectively, vanish. With and Xj^’ both as p orbitals i 

or both as d orbitals or with either as s and as d I 

orbital, this integral does not vanish. The p-p terms are the I 
most important ones, and in an sp basis set, when the d orbitals i 

■ ' I 

are absent, they are the only ones to be accounted for. 


involving s and s orbitals or s and p orbitals as 


Using the well-known forms for the hybrid orbitals, i 
expressions^ * for Townes-Dailey approach! 
are available. The EFG*s are expressed in terms of electron | 
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populations on the p orbitals of the atom in question. All 
these expressions also involve»the field gradient due to the 
p~crbital of the free atom»which is used as a parameter. Fur- 
ther, conceptual parameters such as hybridisation, ionic chara- 
cter and pi character enter the expansions, making the predi- 
ction cf EPG’s by this approach less objective and too empirical. 
The contributions from the inner core of atoms are completely 
neglected. Thus the Sternheimer effects ^ are ignored. Some 
attempts have beep, made to remedy some of the drawbacks in the 
Townes-Dailey approach. Particular mention may be made here of 

Mr 4 n* 

the work of Cotton and Harris , Sichel and Whitehead and 

1 ft 

White and Drago . Although these authors have attosmpted to 
decrease some of the deficiencies of the Townes-Dailey approach, 
they still either neglect some cf the terms in eqn. (1.29) com- 
pletely or approximate them with some rough estimates, and in 
this prcicess heavily parametrise the calculations. It has been 
shown by O'Kcnski and Ha that the approximation used by Cotton 
and Harris for the two centre integrals, is ra-ther unsatisfactory. 
Many authors^’ have employed p-orbital populations 
obtained from semi -empirical as well as ab initio MO calculaticns 
with a view to correlate the EPG's of first and second row atoms, 
especially using the Townes-Dailey type approach and have obtained 

varying success. However, the parameter 1*^)free atom 
in these calculations for a given atom hats varied widely (-7 MHz 
to -14. 1 MHz for and naturally, casts serious doubts on the 
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soundness of the procedure. In fact, very recently, Barber 
20 . 

et al. failed to obtain a correlation between their ab initio 

1 4 

populations and IT EPG's in nitrogen containing compounds. In 
a recent work based on a version of semi-empirical MO formalism, 
referred to as MIlTDO/3, Dewar and coworkers^^ have calculated 
BFG' s at IT-site in some molecules. These authors have employed 
the point-charge approximation for the two-centre integrals and 
the Mulliken type approximation for the three-centre integrals. 
To our knowledge, at the semi-empirical SOP MO level, only two 
groups of workers have made a complete calculation of the 

EPG's based on eqn. (1.29) without approximating any of the 
integrals. However, both these calculations concern themselves 
with EPG's at the deuterium site only. 

1.5 Scope of the Present Work 

Amongst the semi-empirical MO formalisms'^ ,CHD0/2 and INDO 
formalisms have been extensively employed for the study of 
various molecular properties. However, as far as EPG's are 
concerned, the wave functions from these formalisms have, been 
employed mostly in the fram.ework of the Townes-Dailey approach. 
Claims regarding the superiority of these wave functions for BPG 
calculations cannot therefore be accepted without reseivation. 
It is essential that these wave functions are properly tested 
for several nuclei and a variety of molecules, by the evaluation 
of EPG's with the inclusion of all the integrals in eqn. (1.29). 
Such a test has not been carried out so far. Two recent 
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23 24 

calculations ’ have addressed themselves to .this point 
using only one nucleus, namely deuterium. In the present work, 
we hove calculated the BFG' s at DjByN and 0 sites in about 
twenty five molecules using CNDO/2 and INDO wave functions with 
Slater exponents and including all integrals. ¥e have also 
examined Cinother version of the HTDO method due to Nanda and 

p C 

Nara.siinhan which differentiates between orbitals with 

different m-, values such as p , p and p and yet retains 
i- X y z 

rotational invariance. It is well known that these semi- 

empirical formalisms with the neglect of differential overlap 

correspond to calculation over orthogonal atomic orbital (01.0) 

27 

basis sets" . Hence, in order to calculate the EFG' s we have 
assumed that these semi-empirical MO's are over a non-overlapping 
basis and transformed them to an overlapping basis. This pro- 
cess is often referred 'to as ' deorthogonalisation' . Ill our 
calcula.tions have been performed with deorthogcnalised basis 
sets. Our results have been compared with ab initio results and 
experimental results, wherever available, in order to bring out 
strengths and weaikenesses of these three semi-empirical formali- 
sms. The evaluation of BFG's in molecules requires knowledge of 
the ground state wave functions only. Since the EFG operator is 
a one-electron operator and samples the region of space close to 
the nucleus, the behaviour of the semi-empirical wave functions 
in this region can be also properly assessed. The role of or- 
bital exponents in the STO's used can also be examined, by their 
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on the calcsfiJilated IFG-'s, HEhanks to th® rork of Snyde^^ 

28 

and 3ia-sch vfe ha'»^5 a compilat iOQ. of ah initloct results^ on 
EPG" a for several laocolecules, ohtKiined using dO‘Ui.t}lle zeta STO’’® 
in felh® basis set- We are thus d;n a fortunate jittsition to cotti-^ 
pare- uniforiQly the? ®'erf orraance a:f the three se'a.la—einperical KGO 
fornaDisms rnentio-W'd earlier, 3ig;sainst ab iaitS.ro methods in 
most ©ases. 

in Chapter IE af the thesiss wre present an. outline of the^ 
sernt “'empirical fon-ais.lisms emplojr'ed and also dsa’Cciribe in detail 
the procedure for the evaluation of the EP3- laftssgrals. 

2 

Chapter III present's and analys&si our results 'O^iri EFG-'s at H, 

p~nd '^'^0 sfltttes in some liciea.r inolecul6'.Si'Wb.ile Chapter HY 
pres^'eants and airalyseis EPG results in some no 3 H“-.lIinear as well ( 0 iS 
eye Ll«s molecules. (CSlhapter T presents further rfieeults invo Ivii.n® 
D,B,B8 Eind 0 sites .iiiuo several typ icoal molecules., The thesis 
0 on.cliuides vith soms'S remarks end siunggestions foriTuture work. 
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CHilPTER II 


ALL- VALENCE ELECTRON SEMI -EMPIRICAL MO FORMALISMS 
ANL CALCULATIONS OF ELECTRIC FIELD GRADIENTS 
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II. 1 G-eneral 

The present calculations have been performed using- -ttiiree 

different semi-empirical all-valence electron ICA.O-MO fc-raalisms , 

all of which fall into the categcryr C'f what is popular Ijr linown 

as ’Keglect of .Differentia.! Overlap' (NDO) methods^'^, The 

nomenclature de.rives from the fact tliat certain two -electron 

integrals in the Hartree— Fock— Roo tbas,n equations^ are neglected; 

this was paraphraised in the origina.1 literature as' differential 

overlap is neg liscted' in these intesgrals. Hcwever, it has hoen 
5 

shown that what these approximations dc»is to simulate a cal- 
culation using a basis set cf ortho go nalis '^d atomic orbitals 
(OAO). Thf-^ more popular of th'^se methods have been C_on]plete 
Neglect of Dif f erentisil Overlap (CUDO) and Intermediate Neglect 
of Differential Ou-erlap (INDO) schemes'^ . Several schem«es of 

integral evaluation and approximations have teen implemented 

6 

at these levels . For example, a modification of the IffDO 

method known as the MINDO method has been developed by Pevar 

and cowcrkers. Several versions of KINDO are available, WINDO/3 

7 

being the latest Another method tbe MGZDO ( many-c ent r aero 

If 

differential over .Lap) has been d.ev(=lo»ped by Brown and Roby . 

The most sophisticated of the NDO ni«thods is the NDDO method, 
originally developed by Pcple amd covrorkers^^» but computet Lons lly 
implemented much later. These methods take into account only 
the valence shell of electrons, it being assumed that the Inner 
'core' of electrons is una.ffected by bend formation. This 
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assumption, it must be pointed out here, has serioue iHimplica- 
tions in ihe calculation of nuclear quadrupole coupliijtag con- 
stants (JQCC) from theoretical values of EFG- evaluate(il0 using 

Q 

'valence only' schemes. It has been shown by Sterniteiftier » 
that the nuclear quadrupole moment can polarise the' sswirrounding 
spherical charge distribution and lead to field gradi<«;giit contri- 
bution frorat closed -shells around the nucleus in questii&n. This 
effect known as * Sternheimer Shi elding/ Ant is hi elding Effect' ?5an. 
be accounted for in the all-electron calculations. Ho'owever, in 
'valence only' calculations it is generally assumed tfe^nat the 
Sternheiner Effect is either negligible or is constawntt in a 
series of related molecular environments. ¥e shall r^gtturn to 
this q.uestion later in this thesis. 

¥e have chosen for our studies the CFDO/2 and IllDOD methods 

p TU 

of Pople and coworkers'^ , and a method at the INDO lievv&l using 
an explicitly orthogonalised atomic orbital basis d ew^glop ed hy 
landa and Haras imhan^ (M-BTDO). 

II. 2 ICAO~MO-SQF Methods and Atop roximat ions 

Ve restrict ourselves to the closed -shell case, ttihere the 
total wave function is expressible as a single deterraiiaant of 
doubly occupied molecular orbitals, and define our basgiis as the 
valence basis, for such a molecular system, the Hartu'-ee-Tock-. 
Roothaan equations^ are 


( ir.i) 
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where N is the number of basis functions and are the ICAO MO 
coefficients. Here is the overlap integral, and is given 

by 

^ ^ P^j^[(i3/kl) -^(ik/jl)] (II. 2) 

k 1 

where is the one-electron cor e-Hamiltonian matrix element, 

^kl charge density matrix element given by 

occ. 

^kl ^ I ^mk ®ml (II-3) 

m 

and 

(ij/kl) = J {<>^(1) ♦j(l) i .>i^( 2)‘>3^(2) dr^ d?2 (II. 4) 

Normalisation requires that 

I S.j = 1 (II-5) 

i J 

The set of equations (II, 1) are to be solved iteratively, since 
the F matrix elements are themselves functions of the c's via 
eqns. (II.2) and (II. 3). The MO’s thus obtained are popularly 
known as SCF MO's. The major stumbling block in the solution of 
eqn. (II.l) is the number and variety of two -electron integrals 
to be evaluated. However, if the basis set is orthogonal many 
of these integrals become zero or very small and may be neglected. 
It is observed that these integrals which become very small are 
the ones involving an overlap charge distribution of the type 
‘!*i(l) with i j , so that ’neglect of differential overlap', 

viz., neglect of 4^(1) ’^j(l) leads to a drastic reduction in 

the number of integrals. At the CNDO level, all integrals 
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involving an overlap charge distribution are neglected; at the 
lUDO level, one-centre exchange integrals (ij/ij) are retained 
in addition to those integrals retained in the CNDO level; and 
NDDO retains, in addition, all integrals not involving a ’diatomic 
overlap charge distribution. Also, because of the simulation of 
the effect of ah OAO basis>S in eqn. (II. 1) may be replaced by 

-L J 

6 . . >so that the set of equations to be solved becomes 


I 

i 


c . 
mi 



B 6. . ) = 0 
m ij 


(II-6) 


and the normalisation condition becomes 



(II. 7) 


Different approximations and parametrisations are employed in 
different semi-empirical methods. Thus, each semi-empirical 
method is characterised by the way the core-Hamiltonian matrix 
elements . and the remaining two-electron integrals are evalua- 

X J 

ted. The overriding concern in semi-empirical MO theories in 
the choice of the parameters is the so-called ’invariance' pro- 
blem. The exact solution to the Hartree-Pock-Roothaan equations 

4 

has the property that the calculated results remain unchanged 

when orthogonal coordinate transformations are made. However, 

when approximations to the Pock matrix are made and some integrals 

are neglected, this property need no longer be satisfied automa- 

1 2b 

tically, Pople and coworkers / address themselves to the 
problem of 'rotational invariance’ and their prescription to solve 
this problem is to use values for two-centre coulomb integrals and 
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off-diagonal core-Hamiltonian elements dependent only on the 

2i1d 

atoms concerned and not on the individual orbitals . Pople 

suggests that the results be invariant under hybridisation also. 

1 0 

Jug , however, is of the view that invariance requirements can 
be met by employing even less drastic conditions than demanded 

5 , 

by Gi\lD0. Cook et al. emphasise that invariance is to be demanded 
of calculated properties and not of the integrals entering the 
calculation. Cook is also of the view that the widely used 
invariance 'principle* of approximate MO theories is physically 
unreasonable and formally unnecessary. However, approximations ^ 
of coulomb integrals by integrals over spherically symmetric 
S“Orbitals> in order to ensure invariance, as is done in CNDO and 
IKDO, may be a drastic step in that, it neglects important spatial j 
distinctions among the highly directional p- and d- orbitals. 

This is likely to be of serious consequence as far as tensor pro- ; 
perties like EIG are concerned. 

We shall now examine in detail the expressions for the matrix 
elements in the three semi-empricai formalisms chosen by us, and I 
indicate the relevant approximations and parametrisations employed; 
The diagonal elements of the core-Hamiltonian, , may be written 
as 

= (i I - ya I i) + (i I - Vg j i) (II . 8) 

where A is the nucleus on which orbital i is centred. It is seen 
that the first term is essentially 'atomic' whereas the terms 



27 


(i j- Vg ji) are of the nature of ’corrections' for the presence 
of other nuclei. It is, therefore, usual to set 



U. . 

11 


I, V 

Bt^A 


(II. 9) 


C NDO/2 Paranxetrisation 

In the CK'DO/2 method 


the approximations and parametrisa- 


tions used are the following : 

i) All exchange integrals, and three- and four- centre integrals 
are neglected. 

ii) The coulomb integrals are reduced to one per atom pair, and 

hence denoted as They are evaluated using the s-type 

AO's on the respective centres. Thus, no differentiation 
is made between s,p,d, .... orbitals. 

iii) Off-dia,gonal elements of H are taken to be proportional to 
the overlap integrals, the proportionality constant being 
characteristic of the atom pair. 

H. . = (S. j/2) + Pj) (11.10) 

where A,B are atoms on which orbitals i ,3 respectively are 
centred. The p's are chosen empirically to obtain agree- 
ment with ab initio results. It may be pointed out that 
it is here that OAO simulation enters the theory, the off- 

diagonal elements given by eqn. (11.10) are just of the same 

5 

order of magnitude as over an OAO basis . 
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iv) is evaluated as 


U., =- ysd. + A. ) - (Z^ - 1/2) 


( 11 . 11 ) 


where stand for the ionisation potential and ele- 

ctron affinity of the orbital i respectively and are 
obtained from atomic data. is the core charge for 

atom A and is the one-centre electron repulsion 

integral. 

B 

v) is evaluated as 


A® = - "ab 


I lMDO lara m etrisation 


( 11 . 12 ) 


1 3 

In the INDO procedure of Pople, Beveridge and Dobosh , 
as pointed out earlier, the major difference from CNDO is that, 
one-centre exchange integrals are retained. These are evaluated 
using Slater-Condon B and & parameters and G integrals 

are obtained from spectral data while B° value is obtained theo- 
retically as in CNDO, using Slater s-type AO's. Also, because 

2 1 

of the introduction of B and G terms the precise expressions 

for U. . are somewhat different from CNDO. 
ii 

Nim -INDO 1 arametrisation 


The method of Nanda and Narasimhan^ (NN-INDO) is at the 

INDO level, but makes explicit use of lowdin-Orthogonalised ■ 

orbital basis. The orbital exponents are chosen by Burns* 

1 6 

rules . Nanda and Narasimhan consider the transformations of 
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AO’s under coordinate transformations and arrive at the rota- 

Q 1 7 

tional invariance criterion ’ for the one-centre coulomb 
integral J as 


^ii 


= J. 


11 


+ 2 


^ii 


(11.13) 


where i and i’ are orbitals corresponding to the same 1 value 
and , is the one-centre exchange integral. Thus, Nanda and 
larasimhan are able to differentiate orbitals with different 
values in repulsion integrals of their lEDO method, and still 
achieve rotational invariance, whereas in Pople's INDO pro- 
cedure such a situation does not obtain. Also, Barns' rule 
1 6 

exponents employed in MA-INDO are different for s and p 
orbitals, unlike the equal s and p Slater exponents used in the 
CNDO and IHDO methods. 

One-centre integrals and the one-centre two-electron 

integrals) are derived from atomic data"'^'^^ in such a way that 
eqn. (11.15) is satisfied; there is no equalisation of all cou- 
lomb integrals as in the case of CADO or IBDO of Pople. The 
quantities thus obtained are partitioned into a potential 
energy and a kinetic energy component as 

U. . = T. . + lA (11.14) 

IL 1 IL 1 

by making use of T^^^'s obtained by a suitable scaling of theore- 
tical kinetic energy integrals over atomic orbitals orthogonali- 

A 

sed to the core AO's. The 's thus obtained by difference are 
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used for obtaining the framework corrections to the core- 

"R 

Hamiltonian (see eqn. II. 9) using the relation 

A 

\ ~ ^ - (l-exp(-2 ^ (II. 15) 


which is easily seen to have the correct asymptotic behaviour. 


In the above expression 



is 


the average repulsion integral 


between orbitals on A and orbitals on B, calculated using Ohno's 
20 

relation 


/■•s in a.u. 


(II. 16) 




The atomic Y' s are obtained by the method of Yamaguchi and 

21 18 
fueno used in conjunction with O'Leari's atomic data . The 

two-centre coulomb integrals are also obtained from the above 

relation as 


(ii/jj) = (11.17) 

[Ri^^ + (2/(ii/ii) + ( 03 / 33 ))^]^ 

The off-diagonal core-Hamiltonian elements are obtained from the 
relation 


y2(H, 


+ H, .)s, . 

3 3 13 


y2(T, 


"33^^13^^ 


(11.18) 


with the term S. .(1 - IS. -j) evaluated over a diatomic coordinate 

X j ^3 

system and then rotated to the molecular system to maintain rota- 

22 

tional invariance . The Hamiltonian over OAO's is obtained by 
transforming the H-matrix as 


H = S~ H 


(11.19) 
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where the bar over H signifies OAO basis. The two-electron 
integrals from atomic data are used uncorrected, since it is 
known that they are not as severely affected by the OAO trans- 
formations as are the core-Hamiltonian elements . The atomic 
parameters used in the M-INDO calculations for the first row 
elements are given in Table II. 1. The parameters for lithium 
were not readily available to us and hence we obtained them by 
extrapolation using a cubic polynomial fit of the values for 
Be, B, C and N as a function of atomic number. 


I I .5 Methods of Evalua t ion of Electric Field G-radient Integrals 

For the evaluation of the EFG- tensor components using an 
LCAO MO wave function it is necessary to evaluate integrals of 
the type 


( ^ A Ijl I ^ N 


( 11 . 20 ) 


•A 


where r^^ e ( ^ components of the radius 

vector connecting nucleus A (at which the EFG is desired) to a 
given electron, and (j)^ are the AO's used in the calcula- 

tion. This gives rise to the following four types of integrals 


(*1,^ I Oj^l ) 


I.YPe I : (*v“ 10,1 (II.21a) 

with both orbitals on the same centre as the operator (one- 
centre integrals). 

Type II ; (<»/ ( Qj^ 1 »/) (II. 21b) 


(two-centre integrals) 




Atomic Parameters for the Pirst-Row Elements Used in the NF-INDO Method 
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( 11 . 210 ) 

(two-centre integrals) 

^xmzLl - I oj (i)^°) (ii.2id) 

(three-centre integrals). 

Type I integrals are easy to evaluate in closed form over 
the STO basis set by separation of variables. Bdlund et al. 
have listed formulae for some of these. In the past, calcula- 
tions of El'G using the NDO methods have mostly neglected the two 
and three-centro integral contributions to the EBG- operator. 
The importance of two- and three-centre integrals in the evalua- 
tion of EPG- ' s by semi-empirical MO theories was recognised by ■ 
Dewar and coworkers . However, in their MIl\iD0/3 studies 
they used only approximate values for these integrals. In the 
case of EB'G calculations involving deuterium? there is no justifi 
cation to discard the two- and three-centre contributions since 
the Is contribution to EPG is zero by symmetry. Hence, to have 
a consistent scheme of evaluation of EEG using semi -empirical 
EDO methods»we need to incorporate the two- and three-centre 
integrals. To our knowledge there are only two studies^'^’ 
reported in the literature where all the two- and three-centre 
Bi'G- integrals have been evaluated in a semi- empirical MO calcu- 
lation, However, these two studies have concerned themselves 
solely with deuterium. 

The two- and three-centre integrals need careful handling, 
especially with regard to the diagonal components of the BEG 
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tensor. Because of the singularity of the integrand at the 

nucleus, the integrals are conditionally convergent for some 

functions and the value of the integral will depend upon the 

way the region around the singularity is handled^^. Stephen 

and Auffray noted this in their work on HD"*", and evaluated the 

integral hy a numerical integration method with a spherical 

volume around the nucleus excluded; they thus obtained the limit 

as the excluded volume tends to zero. Subsequently it was 
28 29 

shown ' »that the same effect is achieved by adding to the 

integrand a term equal to 4 tc/ 3 times a delta function centred 

at the singular point. This has been noted also in the work of 

30 31 

Kolker and harplus and Kern and Karplus . McConnell and 
32 

Strathdee failed to take account of this in their work on esr 

29 

hyperfine coupling as was pointed out by Pitzer et al. 


Por Type II integrals methods based on prolate spheroidal 

33 

coordinate systems or the Barnett-Coulson expansion may be 
used. In the methods based on the prolate spheroidal coordinate 
systems, these integrals (as also others having r ^ in the inte- 
grand, with n > l) necessitate new auxiliary functions in 

34 

addition to the usual and B^{b) functions . Kolker 

and Karplus^*^ generalised Ormand’s Bu^(x,y) functions^^ to 


define 

B^(a,b) 



1 

J dp exp ( -ax 
1 


-b kJ-) 






m 


( 11 , 22 ) 


which along with the A^(a) and Bj^(b) functions lead to the 
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required integral formulae. Huo and Roothaan^^ developed a 

method based on auxiliary ffunetiohs related to Rodthaan’s 

Y ^ 5 , t 

details are not available in the literature, 
Barfield and co-work- ers^^ have published formulae for integrals 
involving 1s orbital on centre A (vide eqn..II.21b} and Is, 2s 
and 2p orbitals on centre B based on the expansion method- of 
Barnett and Coulson. 

Type III integrals were evaluated by McConnell and 

32 ' j 

Btrathdee who failed to take into account the d-function 

and thus were led to wrong results, as was pointed out by 

29 

Pitzer et al., . Pitzer et al, made use of expansions of the 

•xn I 

operator on centre A onto centre B following Hobson^ . Recently, j 
23a 

Bdlund et al. evaluated these integrals but Barfield and 
25 

coworkers have pointed out some errors in their results, and 
have also listed explicit formulae invol-ving orbitals upto 2p 
based on the earlier work of Barfield^®'^^, 

I'or Type IV (three-centre) integrals the only available I 

method seems to be the ’Gaussian-Transfoim Method* of Sha-vitt 
and Karplus'^*^'^^ For the BFG integrals this method was adopted 

4 1 

by Kern and Karplus"^ . The Gaussian-transform (GT) method has 
the advantage that in the limit as the nuclear centres are made to | 
approach each other, the formulae go over to the correct two- 
centre formulae. We have in the present work evaluated all the | 
integrals other than one-centre integrals by the GT method. In j 
spite of the fact, that specialised methods are computationally 
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cheaper for the two-centre integrals we have employed the GT 
method for these also because of its ease for computer progra- 
mming, In the next section we give an account of the GT method. 

.I.I.j-4, The Gaussian Transform Method for Evaluation of EPG 
Integrals - . . 


Here use is made of the facts that many centre integrals 
involving gaussian functions are easy to evaluate, and that STO's 
can be expressed as integral transforms of gaussian functions^*^'^^. 
For example, the Is-STO may be written as 


e 


er C 


J 3-5/2 ^/4S 



(11.23) 


o 

The method essentially consists in obtaining the expression for 
the integral involving the Is-type GTO’s and then obtaining the 
Is-STO formulae therefrom by applying the transform (eqn. 11.23) 
to each GTO used. The basic Is-STO formulae »then» become the 
starting point for derivation of the higher orbital formulae by 
the differential operator technique. ¥e outline the method in 
three stages as indicated above, viz., 


i) the derivation of Is Gaussian formulae, 

ii) obtaining the Is exponential formulae therefrom and 

iii) generation of higher orbital formulae. 


The treatment given here for the diagonal components of 

31 

BFG follows Kern and Karplus"^ with some notational changes. 

The expressions for the off-diagonal components have been derived 
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by the author along the same lines. For convenience in for- 
mula generation we express the formulae in terms of the 
functions rather than the s used by Kern and Karplus. 

II. 4. 1 G-aussian 'Is' Formulae 


For purposes of illustration »we shall consider in particular 
two of the BFG integrals, namely, those involving the and 
q operators. The problem of singular integrands occurs with 

jQf 

the diagonal components q , q and q . It has been demon- 
strated that the physically correct value is obtained by adding 
a term equal to 4 ti/ 3 times a delta function about the singular 
point for these integrals. Thus we may write for operators 
centred at C 


3z, 


'C 


^ 6(rQ ) 


^zZfC 

and 

‘^xy,C 


9z, 


(11.24) 


'C 


yg 

5 


9 ' 


r^' 


3 y_ a x„ ^ ^ 


(11.25) 


Let us consider the three-centre case, where Is gaussian 
functions on centres A and B, assumed distinct from C, are 


involved. For q 


zz.C 


we get using eqn. (11.24) the relation 


Ig = I dr exp(-a^ r_^ - ag^g) 


3 z; 


(1/rn) + 


— 2 2 
411/3 f dr exp(-a 2 ^ - a2 r^) 


(11.26) 


Because of the distinctness of the centres A, B and C , the 
operator 3^/3 z^ leaves the exponential part in the first term 
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unaffected, so that we can, after interchange of integration 
and differentiation operations, write 
2 

I 5 = -- 3 ' " S'!? eip(-a, 4 - “2 

4tc/ 3 Jdr exp(-a^ " «2 (11.27) 

Now the integral in the first term in the r.h.s. of eqn. (11.27) 

is the three-centre nuclear-attraction integral, which has been 

evaluated through the use of the auxiliary functions^^ 

1 

Fjj^(t) = I exp(-tu^) du 

° (t > 0 , m = 0 , 1 , ...) (11.28) 

which obey the recursion relations 

= -ft (11-29) 

and 

V") = dr (2t . e-h (11.30) 

The second term in the r.h^s. of eqn. (11.27) may also be 
evaluated using the same auxiliary functions. We define 



x« - x.„ with similar relations for a ,a ,• 
A B y z 


o2 , «2 ^ „2 

a + a + a„ ; 

X y z ' 


(aya.+a2) -t- (a2/«i+°f2^ ' 


(11.31) 


( q .j / q! I +0C2 ) + ( (12/ "I ■*■^2 ^ Njg > and 

(oc ff ^ ^ ^ ^ -j 2 ^ ^B 
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and PC as the length of the vector joining (P^,P„,P_) to C. 

X 

This leads to 


c, o OC OTo n 

^Gr ~ ^ ^ ^ 2 ^ ^o^ ^ ^ ®3cp[- AB ] } + 


4tx {P 3 _[(a^+a 2 ) PC^] -(2/3) (a^+a 2 ) PC^ P 2 [(a^+a 2 ) x 


a. ttp p 

exp(- AB ) } 

“2 


(11.32) 


Performing the indicated differentiations using e^n. (11,29) 
leads to 


(«1+a2^ J’2[(ai+a2^ 

a, ap ^p 


(11.33) 


Similarly for the operator we get 

2 

J ex?(- oc ,4 - agi-^) -^Vy- ( l/i-g ) d? 


T^n;- f “l 


8 'n:(a^+a 2 ) PC^ PC^ P 2 [(a^+a 2 ) PC'^] x exp[ - AB"^] 


1 ^k" ^2 


OC 4 OCp p 

] A T)^" 


(11.34) 


Now let us consider the two-centre case. Here let us 


assume that centres B and C coincide. ¥e may write then 

2 

^&,B=C = Jexp(- a:, 4 - cj r^) (1/rg) d? 

2 

= / { lim 8 xp(- 01 , r? - <x~ rh} — ( i/r ) dr (11.35) 
B ♦ C 3 Zq _ 

Assuming that the limit and integration operations may he inter- 
changed, we now get 
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4 - “2 4) t4- 

®c 

So that the two-centre formulae may he obtained as limits of the 
general three-centre formulae. 

II.4 .2 Formulae Over 1s-3TQ's 

As before, we consider the case of the operator. Using 
eqn. (11,23) we can write the product of two Is-functions cen- 
tred on A and B as 




1 -"A ^2 -^B' 

^2 “ ^ , h-3/2 

~Tr f J (“1 “a* 

0 0 

2 2 

exp[-(o:^ r^ + a2 r^)] da^ da2 




1 "2 
(11.37) 


so that the required integral becomes 

3 - r^ 

I - f ^ 0 

' (I C 

I'lZf 


(1s)a (1s)g dr 


c, c 


C 2 2 

i da^ I da 2 (a^ e3cp[- ^(^ + ^ 


a. 


0 


^ , 2 ,2 

^ 3 2 2 

k- exp[ - - ^2 ^B^ ^ 

■ 


(11.38) 


where the factor within the braces is easily seen to be in 
eqn. (11.26) and which has been already considered in detail in 
Section II. 4. 1 . 

Substituting for I(j from eqn. (11.33) and defining the 
following variables 


a 


<*, + Qig I 


u 


i 


«1+“2 ' 
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X 


u X, + (l-u) Xt, - 


B ~ ^ 


aad other similar 


equations for 


Py» Pg f 


w 

f 

g 


P 


+ P,r +■ P 


u(l-u) AB^ 

i fil + _k>i 

4 'u 


(11.39) 


and 


leads to the expression 

1 


"E 


3 


I du j dt [u( 1 -u )] 


t”"^ exp (-ft 


- g/ 1 ) X 


(3P2 ~ P^) I'2(wt) (11.40) 

We note that AB in the above denotes the distance between points 
A and B, 

Using the auxiliary function 

oo 

Jn,m(w,f,g) = I dt t"^"^ ®xp(- ft - g/t) (II..41) 

o 


the expression for 1^, reduces to 


2 C, 


1 


^2 2 

-P 


"E 


4 — ^ f du ^ ^ Ji 5 

5 O [u (l-u)]^^2 


(11.42) 


Now the integral evaluation can be done provided the function 
Ji o os-ii evaluated. This can be done using the known series 

■ffC 

expansion for ^^^^(t) in eqn. (11.40). Another approach to the 
evaluation of the J„ _ functions is via the related functions 

XJl ' f SJU 
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m-lH-4- 


w 




(11.43) 


where 

0 = f/w, X = gw 

¥e have 




r(m+4' 


y j P^(t) exp(-0t - dt (11.44) 


Substituting for Fjjj(t) from eqn. (11.28) and reversing the 


order of integrations leads to 


t 


T, Ja,x) = — ^ 

r(mH)(2'i:) 


° (11.45) 


where k^(x) = K^(x), K^(x) being the modified Bessel fun- 

42 

ction of the second kind . The integrand being a smooth fun- 
ction, evaluated by gaussian quadrature. In 

the present work the functions K^(x) are evaluated directly 

42 

for n = 0,1 and through the recursion formula 

+ ^ V-) 

for higher values of n. Kq(x) and K.^(x) are evaluated from 

42 

their series expressions for 0 < x < 2 and by a polynomial 
approximation^^ for x ^2. The reduced functions k^(x) are 
obtained from the K^(x) and used in the evaluation of the 
and thereby "the 


functions needed. 
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As for the EFG operator, we have similar ly 

1^, = 2 ?2 i 


r — Ji «(w,f,g) 


(11.47) 


with the same notation as defined in eqn, (11.59) 


II. 4. 3 


er Exponential Orbitals 


Relations such as 


"■a ® 


X'l^ © 




r r 

■2 B 


3 " ^A 

3 1 a 


^ ^2 ^2 


^2 




etc. (11.48) 


are used in generating the formulae involving higher exponential 

orbitals. One simply differentiates the Is - formula with respect 

to the parameters appropriate in each case. However, this tends 

to be a tedious and error-prone task if done manually. Hence, 

these formulae have been generated by us automatically using a 

computer program with the aid of the ’replacement operator 

40 41 

technique'. To illustrate this technique ' , consider a 

general term in the formula 

^2 1 -3/2 + n. -3/2 + n^ i ^2 ^3 

0 = ?:2 i ^ P Px Py Pz ^ 

0 ■ j 

m. m, mT- ' ^ I 

a ^ a a J (w, f,g) du (11,49) | 

X y z n,m^ ' 

with the notations as defined in eqn. (11.39). ^ 

3 ' 

The effect of an operator such as - r~ may be found using the ; 

^ : I 

chain rule, which generates a number of terms with modified j 

values of the indices a^,a 2 , ...» etc. That is, we may look I 
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upon this operator as a sum of ’replacement operators' which 
modify some of the indices, a.nd multiply the resulting terms 


with appropriate numerical coefficients, 


’9 © 


33.4 —1 




Here we have used the relations 


Por example, 
n+1 


we have 
(11.50) 


(w,f,g) 


3 w 




3 <T ( W , f , g ) 

3? 




(w,f,g) 

n.m 

•"i-g 


'^n+1 


(11,51) 


which follow from the definitions of the J 's and eqn. (11.29). 

33. in 

A given higher orbital formula is obtained by successively 
applying appropriate replacement operators. This can be coded 
for a digital computer representing each term by a coefficient 
and a set of numerical values forming components of a vector, • 
like terms are collected out of the resulting terms and the 
appropriate coefficient and indices printed out. The formulae 
so generated are coded for computation, 

Por the purpose of formula generation the formula is 
represented by a sum of two terms, as 

( 1 S l<J^^ I 1 S) = I, + Tj 
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where ^ ^ 

! [u(l-u»5/2 

and ° (11.52) 

1 2 

O L V yj 

In this case, the terms generated from both and T 2 

added to give the final result. For the q EFG- operator, the 

xy 

Is ~ formulae consist of only one term and hence the formula 
generation is more straightforward. Formulae for all the com- 
ponents of the EPS operator have been similarly generated. 

II . 5 Testing of the Computer Program for Svaluation of Integrals 

All the necessary integrals for the evaluation of SFG's have 
been obtained in the present work by means of computer programs i 
written by the author himself, in FORTRAN 10 and have been success-; 
fully executed in a DEC 1090 computer at the Indian Institute of 
Technology, Kanpur, India. Double precision has been used throughout 

Two versions of the programs were created, the first using 

16~point gaussian quadrature both for evaluation of functions j 

and for the final integration, and the second using 12-point I 

44 

gaussian quadrature instead. Since, Scrocco has reported all 1 
the BPG integrals for the case of CO with orbital exponents from | 
Slater's rules, we used the same for checking our results. The j 
integral values from the 1 6-point are compared with his results j 
in Table 11,2. The agreement is perfect to the last digit { 
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reported by him, in all but the cases of (2p^Q | q| 2p^^), 

^^^00 I ^ I 1 ^ where there is agreement 

within 10 o.u. In Table 11,3 we give the integral values 


obtained with the 12-point formulae. It can be seen that the 
values given in Tables II. 2 and II. 3 agree to within 10“^ a.u. 

We have therefore employed the 12-point version in our calcula- 
tions since this results in a near-halving of computer time. 

25 

Barfield et al, have published values of the integrals for the 
components of the BPG tensor in the case of HCU with the expo- 


nents Z 


H 


1.24, Cq = 1.72, = 1.95 and distances rQg.=1.063 A , 


r„., = 1.5553 A, We calculated q and q , for HCN using the same 

wiM zz 

exponent and geometry as employed by Barfield et al. ¥e used the 
12 -point v-^rsion of our program. The results are compared in 
Tables 11,4 and 11,5. It can be seen from these tables that the 
agreement is excellent. We have also tested some two-centre 
integrals using the method of Pitzer et al. ^ and the method of 
Ko Ik er and Karplus^^. 


A listing of the EFG integral evaluation program written by 

the author is given in Appendix I. A table of conversion factors 

between SPC values in a.u. and nuclear quadrupole coupling con- 

2 7 

stanta (NQCC) values in frequency units for the nuclei D, Pi, 

"'^B, and is given in Appendix II. 

In the present work, we have calculated the EPG’s 
for a variety of molecules using the semi-empirical MO formalisms, 
namely, croO/2, INDO and OT-IITDO. Our results on molecules of 
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2 7 11 

interest containing the quadrupolar nuclei D, Li, B, 

14. 1 **7 

nnd 0 fire reported and discussed in the next three 
chaptern of th,is thesis. Our calculated EBG values are quoted 
upto five decimal places and the asymmetry parameter values 
up to three decimal places uniformly in this thesis. 
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Ii r.l i.;(»noral 

In the case of linear molecules the calculations of EPG-’s 
"becomes a much simpler proposition than in the non-linear case, 
since now the molecular symmetry axis and the principal axis 
of the F'lFG tensor coincide. Iherefore with the z-axis chosen 
alon,^ the molecular axis, the BFG tensor is diagonal, and by 
aymmetry and the Laplace equation we have q = q = - t 1 » 

i.e., the asymmetry parameter p = 0 and q^^ determines the EB'G 
tensor compl'^tely , In order to test the three semi-empirical 
MO forma 1 tame f in terms of their ability to predict the EFG’s 
at various quadrupolar nuclear sites in linear molecules^we 
have chosen the following seven molecules viz. LiD, DP, DCET, 

DCCD, NCCN, and CO. The calculated EFG’s for "^Li, ^"^1, 

^ ^0 in these molecules have been compared with other theoretical 
and available experimental results. The experimental results 
are generally reported in the literature as quadrupole coupling 
constants in frequency units and these have been converted into 
BFG’o in atomic units by us»using the conversion factors given in 
Appendix II, Our calculations have all been performed using the 
experimental bond distances for the molecules under study. The 
values of bond lengths used were taken from literature, and 
wherever possible from microwave determinations. The bond 
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distance data used in our calculations are presented in 
Table Li 1,1. 

Since our calculations make use of a valence basis set 
only for the electronic part, it is necessary that the nuclear 
contribution be calculated using a screened nuclear charge in 
place of the atomic number of the nucleus in question. For 
atoms i'roiTi id to F, therefore, we calculate the nuclear contri- 
bution using the core charge 

Zf = - 2 (III.1 ) 

whoro A| is the atomic number. This procedure has been followed, 

1 2 

for oxampl('!, by Bloor and Maksic . Barfield et al. , however, 
componsut',' for the core electrons, differently. They treat them 
as doubly occupying an uncharged is orbital and add their con- 
tribution to that of the valence MO’s; they then use the A^ dire- 
ctly in the nuclear term. However, at typical internuclear dis- 
tances and using typical Is - exponent values, we find that this 
procedure docs not lead to signif icantly different core contri- 
butions than the point charge approximation used by us. At the 
oxygon centre in CO, for example, this results in a difference of 
about 0.003 a.u, in the final computed field gradient. The point 
charge model, we feel, is more in consonance with the spirit 
of the all-valence electron formalism used. We may note here 
that tho HB'G's reported in this thesis are all in atomic units, 
unless ntherwlSQ stated. 
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Table III.1 


Bond duitances (R in 1 ) used for the Linear Molecules 

under study* 


Molecule 

Bond 

R 

Molecule 

Bond 

R 

LiD 

Li-D 

1.5939 

«2 

IM-I 

1 .094 

DP 

D-P 

0.917 

NCCN 

N-C 

1 .150 

DON 

D-C 

1 .066 


0 

1 

o 

■ 1 .380 


G-N 

1 .060 

CO 

o 

1 

o 

1 .128 

DCCD 

O 

1 

1 ,060 





C-G 

1 . 204 





'‘■L.R. Sutton, "Tables of Interatomic Distances and Configura-_ 
tions in Molecules and Ions", The Chemical Society, London( 1 965) . 
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In ordor tc; study tiiG effect of inclusion of three— 
centre' inttyjruls on the computed field gradient in the approxi- 
mate NIO 1. '.'rtdallr.mn , we calculated the field gradients at the 
deuterium centre in DON and DOCD and at the nitrogen centre in 
HGN using the three models chosen, namely CNDO/2, IFDO and 
KJ)J-IND0, arbitrarily setting all the three-centre Integrals to 
zero. Our results from this calculation are compared with the 
r^unrlts .-'btuined with all three-centre integrals included in 
Tab To I 1 1. 2, Kern and Karplus^ in their calculations using 
ab-initio wavi,’ functions have used the same procedure to eva- 
luate thi' L::i|’,''rtanc'' of three-centre integrals. It is seen that 
tho irffoct is much more pronounced in the case of D than in the 
case 'if N; this is to be expected since in the nitrogen case 
one-contre' contributions are the largest, whereas in the deu- 
terium case the one-centre contributions are zero because of 
tho syminotry of the is orbital. In the total field gradient, 
this effect will be magnified still further, since this is 
computed as a difference of the nuclear and electronic contri- 
butions which are of opposite sign and comparable in magnitude. 
Thus, w.' conclude that for deuterium the three centre contri- 
butiC'ns arc quite significant, whereas for the first-row- 
el onionta these are not so important. These conclusions are also 

3 

supx^arted by the earlier findings of Kern and Karplus . 

For tho sake of completeness, we have performed computa- 
tions rspo.rted in this thesis with the inclusion of all 
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integrals. Since we interpret the NDO wave functions as being 
over an OAO basis set, the MO coefficient matrix has to be trans- 
foruicd Si! O.S t'l give the coefficient matrix over the overlapping 
AO bu.Bis set. The transformed coefficient matrix C_ is obtained 
by the transformation 

0. = C ’ (III. 2) 

In this equation C_' has row indices labelling MO's and column 
indices Is.'bolllng OaO's. If the reverse convention is used, 
then oqn. ([[1.2) should be replaced by 

d = 0,’ (III. 3) 


This procbso is referred to as ’ deorthogonalisation* . All our 
culculatif'ns are done using the deorthogcnalised MO coefficient 
matrix. 

In order to analyse the contribution to the EPG' s atomwise. 


we my write the expressions of the electronic contribution to 
the EFG as^ 

n M M M n 


2 i, ,1, hi <^ki VI 


1=1 k =1 1=1 


M 


i; ““ •!>’ ’ 


(III. 4) 


where k,l are the Indices over the deorthogonalised basis set j 
AO’s (M in number) and i is the index over the occupied mole- | 
cular f;rbitals (n in mmber). The prime over the second suiDmatio3| 
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in "tiie r.h.*s. implies that terms with k = 1 are to he omitted. 
Row on the hand side, the quantity within braces may he 

regarded as ,'.'ross population weighted hy the field gradient 
operstor, qQp,oi "the k ^ AO, i.e. the k"^^ AO contribution to 
the BFG. The overlap contribution is divided equally between 
the orbitals involved, so that the sum over all orbitals equals 
the total electric field gradient. Atomic orbital-wise contri- 
butions reported in this chapter have been calculated by this 
procedure. We do not present the electronic contributions MO- 
wise since there is no uniqueness about the MO's. It is well- 
known that MO's may be subjected to an orthogonal transformation 
leading to differing contributions orbital-wise. Idirther, large 
amounts of cancellations in the contributions to EPG’s from 
various MO's do occur making the interpretation very difficult. 

In the following sections we present and analyse the results of 

our calculations on EPG's in terms of the nuclear sites, viz., 

2 7 1 4 1 7 

H, ^Li, N and 0 in the molecules mentioned earlier, 

III. 2 BFG'o at Deuterium Site in liD, PF, BCR and DCCD 

Table TII.3 summarises our results on EPG's in the molecu- 
les LiD, DP, DON and DCCD at the deuterium sites. ¥e shall dis- 
cuss these results now. 

LiD 

The total IFG values obtained using CRDO/2, IRDO, and 
NR-INDO methods are all seen to be higher than the experimental 
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Table III. 3 



'u n.t the 

Deuterium 

Site in LiD, 

DF, DON, and 

DCCD 

Mo1C:,'Gu1G 


GNDO/2 

INDO 

NN-INDO 

Expt , 


N'uc . 

0.07319 

0.07319 

0.07319 


LiD 

Elec . 

-0.01050 

-0.00967 

-0.00272 

0.0491 


Tot. 

0.06269 

0.06352 

0.07047 


Due . 

2.69033 

2.69033 

2.69033 


DF 

Elen . 

-2.19764 

-2.23234 

-2.05927 

0.5270^^' 


Tot . 

0.49269 

0.45799 

0.63066 


hue; . 

1.13280 

1 .13280 

1 .13280 

0.3009^°^ 


Eli.'C, 

-0. 66B62 

-0.65207 

-0.56599 


Tot . 

0.46418 

0.48073 

0.56681 



Euc , 

1 .10553 

1.10553 

1.10553 


DOCD 

Elec . 

-0.63291 

-0.62369 

-0.53022 

0.2976*''^^ 


Tot . 

0.47262 

0.48184 

0.57531 

(a) 1)00 

lief. 4 





(b) Dee 

Hef. 21 





(c) 

Rof. 26 





(d) Coo 

U.;f. 27 & 

28. 
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value of 0.0491 of Wharton, Gold and Klemperer^ using the 

TOO 1 cK.u ltJ.r beam inugnetic resonance method. However, they are q.uite 

close ic. the theoretical values of 0.0574, 0.0694 and 0.0685 

obtaineil by Iiolkor and Karplus^, from the three different 

levels of ab initio SOP MO wave functions of Ransil^. lOahalas 

and Ncsbet'^ reported a value of 0.0549 using STO's and elliptic 

orbitals# Browne and Matsen^ using valence bond functions and 

an int ’rnucloar distance of 1,6 1 obtained a value of 0.0500 

after avi'raglrg ever the mcleculsr vibrational motion. Bender and 
q 

Davldecn' using a 45-term wave function obtained a value of 

0#0466* ThtJ C’ re mcdGl calculations of Avgeropoulos and Ebbing'^^ 

gave vnlucB of 0.0430 without vibrational averaging. Arrighini 
1 1 

et al. whv) performed a SOP MO with Cl at the Hartree-Pock level 

1 2 

obtaln<;d value of 0.0595. Ditchfield et al. using gaussian 

fjrbitaln report a range of values from 0.0683 tc 0.0695 depending 

13 

on tho typo f'>f wave function used. Docken and Hinze using 

the mu lti-c-.‘n figuration OOP (ICSCP) method obtained values of 

0.0493 and 0.0488 respectively before and after vibrational 

averaging, c'vt-'r the ground vibrational state. At the semi- 

2 

crnpiricuil level, the calculations of Barfield et al. yielded a 
value,' ( ,f 0,0595 by the INDO method as against our value of 
0.0635 by the same method; the difference is due mainly to the 
fact that they treat the ls core on Li differently as discussed 
eorlior. These authors using the optimised STO exponents of 
Hehre, Btowart and Pop le^^ obtain a value of 0.0522. How-r.ver, 
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their use of the original lUDO pa^rameters together with modi— 
fiod exponents is not entirely consistent. On the other hand, 
their study coni irms the importance of the role that orhital 
exponents play in EP(3- calculations. In our studies, the OT-INDO 
formalism using Biirns exponents (which are lower than Slater 
■values) is seen to under- estimate the numerical value of the 
electronic contribution to the field gradient, thereby leading 
to a larger value of EPGr in comparison with the values obtained 
from the (IPDO/J? and TNDO methods. It will be seen that the re- 
sults of these latter two methods, both of which employ Slater 
exponents, are olose to each other. The use of optimised ex- 
ponents of Jiohre, Stewart and Pople (generally higher than 
Slater values) leads to a value of BPG lower than that obtained 
using the Slater exponents. The correlation seems to be unmis- 
takable in this case. 

Our values using the M-INDO wave function (see Table III. 3) 

5 

are close to the ab initio SOP MO results of Kolker and Karplus 
viz. 0,6278, using Ransil's MO function which employs Slater's 
exponents. The CNDO/2 and INDO values are much lower than that 
of NN-IMDO as was found earlier in the case of LiD. Blocr and 
Makslc^ using the CNDO/2 wave function after deorthogcnalisation 
obtained a value of 0.519 compared to our value of 0.4927, the 
difference being due to the different inter-nuclear distance value 
Bloor and Maksic also report a value of 0.580 from their self- 
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conoiDtf^nt charge (SCO) model calculation. The INDO value 

obtained in the present work agrees with that of Barfield 

et alt^ who report a value of 0.4584. At the ab initio level 

several calculations are available^^ ’ ^ These theoretical 

results range from 0.54 to 0.757 and employ different types of 

ab initio wave functions. On the experimental side, Muenter and 
? 1 

Klemperer' have obtained a value of 0.5270 from their molecular 
beam electric resonance study. It is thus seen that the EPG 
results obtained from the three semi-empirical MO formalisms 
are all rather satisfactory, 

DON 

Thf! calculated values at the Groo/2 and INDO levels (see 

Table II T. 3) are comparable with the results of Zern and 

Karplus'^, namely, 0.472, using Slater exponents. McLean and 

22 

Ypshimine's unpublished result quoted by O'Zonski and Ha is 

0.317. O'Konski and Ha using Gaussian Lobe orbitals (GIO) 

obtained a value of 0.365. Harrison ^ has also employed GIO's 

to calculate the EPG's and has obtained a value of 0.3827. 

Ditchfiold et al^^ using gaussian functions obtained values ran- 

24 25 

gim from 0.472 to 0.480, Snyder and Basch » ^ using a wave 
function of double zeta quality obtain a value of 0.3687. At the 
serai-empirical level Bloor and Maksic report a value of 0.481 
using OMDO/2 whereas Barfield and coworkers using the INDO 
method report a value of 0.4792 with Slater exponents. Our 
CNDO/2 and IKDO results substantially agree with the values of ' 
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Bloor &nd M&ksic eind B&rfield et 8.1, Use of Heh.re*“S"te'w8r’t— Popl© 
orbital e xp'^nents In INDO -has been shown by Barfield et al- to 
yiold 8 v«ilut5 ol. 0,410, The best available experimental value^^ 
appears to be 0.3009, Most of the ab initio as well as semi- 
empirical results are seen to be higher than the experimental 
value. The NN-IHUO value, in particular, is higher than the 
CNDO/2 and IMDO values. In the light of our earlier comments 
in the LiD case, we may conclude that Burns exponents do not 
appear to bo optimal for EIPG calculations. 

DGOD 


Our culeulntions by the CNDO/2 and lUDO methods (see 

Table 111,3) give values comparable to the ab initio 'Best Atom 

Molecular Orbital' (BAMO) result of Kern and Karplus^, viz., 

0.473. The NN-INDO result is higher. Harrison^^ employed 

12 

GLO'a and obtained a value of 0.3899. Ditchfield et al. using 


gauasian orbitals report values ranging from 0.473 to 0.483 

1 s 

depending ^n the wave function employed. Snyder and Basch 
using a wave function of double zeta quality report a valup of 


0.3854. The experimental deuterium quadrupole coupling con- 


stant 


Plnkort<-;n"^ ' from molecular beam study leads to 


value rf 0,2976. This value is not very different from that 

28 

reported by Millet and Dailey who used, the liquid crystal 
NMR method. Blocr and Maksic's^ CUDO/2 value of 0.473 and 
Barfield’s^ INDO value of 0.485 using Slater exponents are in 
excellent agreement with ours. With the use of Hehre-Stewart- 
Bople exponents, Barfield reports a lower value of 0.447 in the 
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INDO fcrmalistn. As in the earlier cases the M-INDO formalism 
leads to a. hi^^hor BFG value at the deuterium site. This 
yppi'.'iru t' tu' a ci'nsnnuence of the use of Burns exponents. 

Nex_t.~Ieighbour Contribution to BPG's in BCM' and DCCI) 

vie shall now discuss the atomwise contribution to EPG in 

DON and IXIGD and analyse the importance cf next— neighbour 

ccntribut 1-. n at the deuterium site^ The atomic crbitalwise 

cr ntributi' n t,< th*-' BPG's in jDCN and DCCD at the D site are 

pr'-smtad in Tables Til, 4 and IIIi5 respectively. To see 

cl'’ar1y, th.> importance of next-neighbour contribution, the 

contribution from I) and the neighbouring C alone are summed 

and prt'S''ritod along with the total BF&'s in Table III. 6. The 

total EFO valuGS and the partial sum values are easily seen to 

be close to each other. These trends parallel the findings of 

3 20 

Kern and Karplus' , and Snyder . However, the numbers presented 
in Table I IT. 6 should be viewed qualitatively rather than 
quant itativoly, since the crbitalwise analysis on which they 
arc based Is s 'mewhot arbitrary in the partitioning of the 
overlap contribution. In a few of the tabulated instances it 
is seen that the partial sum of the self- and next-neighbour 
contributions marginally exceeds, the total. This is, however, 
compensated by contributions of opposite sign from the remaining 
terms in the total BPG, It is clear that for the deuterium site, 
the main contribution to EPG arises from the near neighbour atom. 
The near constancy in the BP G value at the deuterium site in 
DON and DOCI) can be thus understood. 
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Table III . 4 

Atomic Orbitalwise Contribution to EFG in DON at D Site 


Gontribu- 
tions from 

CNDO/2 

INDO 

HN-INDO 

if 

Deuterium : 

1 s 

-0.06228 

-0.06030 

-0.03604 

Carbon : 

2 o 

-0.1 5850 

-0*16087 

-0.1 5089 


- 0,15216 

-0*13972 

-0.09260 


%£\ 

# 

0 

1 

-0.15696 

- 0.12473 

Tot , ^ 

-0.45483 

-0.43755 

-0.36822 

*^000. 

0.97860 

0.97860 

0.97860 

Mitrogen 

2s 

- 0,04025 

- 0.04214 

-0.04190 

2 p 

-0.06446 

-0.06520 

-0.06310 

2p 

-0,04678 

-0.04687 

-0.0621 3 

Xot. ‘IqXqq, 

-0.15149 

-0.15421 

-0.16714 

q 

‘nuc . 

0.15420 

0.15420 

0,1 5420 

^total 

0.46418 

0.48073 

0,56681 


^ The nuclear contribution to EFG from Deuterium is zero. 
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Table III . 5 


Atomic Orbitalwise 

Contribution to 

Bi'G in DC CD 

at D Site* 

Contribu- 
tions froffi 

CDDO/2 

INDO 

m-INDO 

Deuterium 1^: 




1 s 

-0.06228 

-0.06136 

- 0.04320 

Carbon 1 i 




2b 

-0,16065 

-0.1 6489 

-0.15137 

2p 

“0,14938 

-0.14958 

- 0.10470 

Pp 

-0.14881 

- 0.13790 

- 0.12290 


-0.45884 

- 0.45237 

-0.37897 

^^nuo . 

0.99530 

0.99530 

0.99530 

Carbon 2 : 




28 

-0.02807 

-0,02840 

- 0.03254 

2p 

-0,05098 

- 0 . 0509 s 

0.04644 

2p 

-0.02478 

- 0,02306 

- 0.02332 

'^'^^*%lec. 

-0.10383 

-0.10238 

-0.10016 

“^nuc . 

0.1021 5 

0.1021 5 

0.10215 

Deuterium 2 ; 




1 s 

-0; 00747 

-0.00753 

-0.00635 

q 

0.00807 

0.00807 

0.00807 

^nuc. 




'^total 

0.47262 

0.48184 

0.57531 


values shown are orbital contributions at the deuterium 
labeiiod 1)^ from the two carbons (0. & C 2 ) and the other 
deterium 

§ The nuclear contribution to BPG from Deuterium 1 is zero. 
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III.l KFG at. Lithium Site In 7AT) 

Our values tor the BFG at the Li site by the CNDO/2, INDO 

and ra-TNDO method are -0.04768, -0.04558 and -0.03308. Kolker 

and Karplus' report a value of -0.04012 using the BAMO wa,ve 

function of Ransil . Browne and Matsen° calculated a value of 

—0.0346 from their ab initio wave function while Cheng and 
29 

Byers-Brown " , who extended the work of Browne and Matsen, 
obtained a similar result. Ditchfield et aL^^ using gaussian 
orbltu.ln obtained values ranging from -0.0206 to -0.0259. 
Arrighini, Guidotti and Tomasi^^ using eleven STO’s and confi- 
guration interaction obtained a value of -0.0429. S-Green^*^, 

who made a Cl culf'ulation with 200 configurations obtained a 

1 5 

value of -0.0429. Docken and Hinze using multiconfiguration 
SCF method obtained a value of -0.0399 without vibrational 
correction. 

The experimental value of Wharton et alf using molecular beam 
resonance technique* is 0.0402. It may be noted here that ours 
are the first calculations on '^Li EFG at the semi-empirical level, 

7 

which include all integrals. In comparing our results on Li 
using the tiiree dilTerent semi-empirical MO formalisms, we see 
that the NN-INDO method leads to a lower value. The reason for 
this Is, that the one-centre contribution to BPG is underestimated 
due to the use of Burns exponents. 
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The ilFCr values obtained in the present work for the nitrogen 
centre in and NCCN are presented in Table III. 7. For 

and TTCN,the orbitalwise contributions to the EPG are also re- 
ported in Tables III. 8 and III.9 respectively. 



In N 2 , Richardson's calculations^ using Ransil's wave 

£ 

function suggested that the electronic contribution to EFG 
from the other nitrogen is nearly cancelled by the nuclear con- 
tribution. This, however, is not observed in the present cal- 
culations (see Table III. 8), in any of the three models used. 

The widest discrepancy occurs in the case of M-INDO. Partial 
cancellation of contributions from the nitrogen p^ and p^ orbitals 
lefado to a smaller overall contribution at that nitrogen centre. 
Richardson obtained values of -0.310 and -1.341 respectively 
using Ransil's BAMO and BIMO wave functions. The sensitivity 

of BFG values to the choice of wave function is marked in this 

22 

case. The result of O'Konski and Ha using gaussian lobe 
basis is -1.086. Dunning^^, using the SOP BIO wave functions with 
various basis contractions giving the best energy, calculated, a 
value of - 1 . 3574 . Ditchfield et alJ^ obtained values varying 
from -0.137 to -0.319 depending upon the nature of the gaussian 
wave function used. The experimental value available from NQl 
studies^^ on solid N 2 corresponds to a value of -1.0254. 
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Table III. 7 


WG ' 

' S 

Nitrogon Centre 

for Ng, HGN 

and NCCN 


Molecule 


CNDO/2 

IlIDO 

NE-INDO 

Expt . 


Blec 

-1 .36501 

-1 .27525 

-1 .54068 



Nuc , 

Tot . 

1 .13171 

-0.23330 

1 .13171 

-0.14354 

1 .13171 

-0.40893 

-1 .0254^^^ 


Bloc 

-1.04301 

-0.92510 

-1 .14919 


HCN 

Hu 0 . 

1 .02615 

1 .02615 

1 .0261 5 

1 

■ 

o 

o 

cr* 


Tot. 

-0,01 686 

0.10105 

-0.12304 


I'llec 

-0.97046 

-0.87565 

-1 .28242 


NCGIM 

rtluc . 

0.88237 

0.88237 

0.88237 

-0.941 6^°^ 


Tot, 

-0.08809 

0.00672 

-0.40005 


(a) 

8oe 

Kef. 

33 

(b) 

8ee 

Ref . 

37 

(c) 

See 

Ref. 

38 
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Table III. 8 

Atomic Urbi trilwise Contributions to El'S in at E Site 


Contribu- 
tions from 

CNEO/2 

INDO 

NN-INDO 

Nitrogen 1 ; 




2s 

0.02765 

0.04889 

-0.03827 

2p 

1 .57524 

1 .57524 

0.90356 


-2.12688 

-2.04692 

-1 .80101 

Tot.q - 

^elec . 

-0.53399 

-0,42279 

-0.93572 

^nuc . 

o 

• 

o 

0.0 

0.0 

Nitrogen 2 : 




2 s 

-0.23438 

-0.241 28 

-0.15871 

2p 

-0.14584 

-0.14584 

-0.08534 

2p 

-0.46100 

-0.46533 

-0.36090 

Tot. q . 

‘elec. 

-0.84122 

-0.85245 

-0.60495 

q 

1 .13171 

1 .13171 

1 .13171 

to tc 1 

-0.23330 

o 

« 

1 

1 

-0.40896 


77 


Table III. 9 

Atomic Orbit alwise Contributions to EPG in HCI at N Site 


Contri- 
butions from 

CNDO/2 

INDO 

NN-INDO 

Nitrogen : 

2s 

0.03730 

0.04940 

0.02560 

2p 

1 .62468 

1.65350 

1 .03424 

2p 

-2.09931 

-2.03953 

-1 .88288 


-0.43733 

-0.33663 

-0.82304 

^nuc . 

0.0 

0.0 

0.0 

Carbon : 

2s 

-0.24293 

-0.25073 

-0.16177 

2p 

-0.05476 

-0.05240 

-0.00924 

2p 

-0.28302 

-0.26069 

-0.08863 

Uleo. 

-0.58071 

-0.56382 

-0.25964 

'^nuc . 

-0.99530 

0.99530 

0.99530 

Hydrogen ; 

1 s 

-0.02496 

-0.02464 

-0.01 531 


0.03084 

0.03084 

0.03084 

total 

-O.G 1686 

0. 10105 

-0.12304 
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A value of -1 .223 has been quoted by O'Konski and and 

presumably this corresponds to the NQR data on solid N 2 
at a different temperature. In comparing all the available 
ab initio values it is found that use of BIMO wave function 
of Ransil and the Dunning wave function yield results closest 
to experiment. This observation signifies the importance of 
properly chosen orbital exponents in calculating EPG-'s. At 
the semi-empirical level the OT-IiTDO value of -0,4005 is seen 
to be a significant improvement on the CFDO/2 and IRDO values. 

noN 

3 

In HON, the value obtained by Kern and Karplus is 

34 

-0,481, McLean and Yoshimine reported a value of -1.204 

22 

from their ab initio calculations. O'Konski and Ha obtained 

a value of -0.854 with the extended GIO basis set. Bonaccorsi, 

Scrocco and Tomasi^^, using 54 STO's reported a value of 

-1.1948, Litchfield et al' quote values ranging from -0.322 

to -0.481 from their ab initio calculations using gauss ian 

36 

orbitals, Very recently Barber et al obtained a value of 

-1.14313 using Dunning's sp basis set. To our knowledge at 

the semi-empirical level ours is the first calculation of EPG 

at site in HCN including all integrals. The croo/2 value 

(see Table III. 7) is extremely low while the IKDO value is 

37 

positive and hence of wrong sign. The experimental value 
is -1.010. The m-IKDO procedure is again seen to improve 
upon the performance of C EDO and IKDO. A closer analysis of 
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the cause for the lower values obtained by the Croo/2 and DJDO 
methods will be now made* Table III. 9 analyses the situation 
in terms of the atomic orbital contributions to the EPG at the 
nitrogen site in this molecule. It will be seen that in CUDO/2, 
the total contribution from the nitrogen centre is nearly can- 
celled by that from the carbon centre. In IlIDO, the contribu- 
tion from the carbon centre is even larger than that from the 
nitrogen centre itself and thereby yields a positiv e EFG value. 
On the other hand* the OT-INDO procedure yields atomic orbital 
contributions which are more reasonable. 

MOON 

In the case of cyanogen molecule* our results (see Table 
III, 7) again reflect the same trend as in HCN. CNDO/2 yields 
a small negative value; IHDO gives an even smaller hut positive 
value; NN-INDO gives a negative value about one-third the 
ab initio result of Bonaccorsi, Scrocco and Tomasi" . These 
authors used a SCP MO with 64 STO's in their ab initio calcula- 
tions. ?ery recently, Barber et al. made an ab initio calcula- 
tion using the Dunning sp basis set and reported a value of 
-1.13716. The solid state NQR experimental data of Haigh and 
Guibo^® leads to a value of -0.9416. Amongst the three semi- 
empirical methods employed here, the ro-IKDO formalism appears 
to be more satisfactory for the SPG at the nitrogen site. 
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Before eoncluding this discussion on the EFG-'s at the 

nitro^^en sites in these molecules , it ma.y be pointed out here, 

39 

that Dewar et al. have calculated the EPG’s by their MINDO/3 
method for a number of nitrogen containing molecules including 
N 2 » HON and NOON. Unfortunately, these authors did not evaluate 
the two- and three-centre integrals exactly and hence no meaning- 
ful comparison is possible between our results and theirs. 

Ill, 5 Field Gradient in CO 

The atomic orbitalwise analysis of the EFG at *'"^0 site in | 
this molecule is presented in Table III. 10 along with the total 
EPG values. 

3 1 

Richardson has obtained a value of -0.151 from an ab ; 

initio calculation. However, an approximate procedure was employed 

40 

by him in the evaluation of the integrals. later, Scrocco : 

6 

reported a value of -O.I 6 using Ransil's wave functions while 
Nesbet^^ , who used a wave function of double zeta quality with 

42 i 

polarisation, obtained a value of -0.627. Lef ebre-Brion et al* 
and Huo^^ obtained with SOP MO procedures, values in the range of | 
-0.621 to - 0.679. Mclean and Yoshimine' , using Slater fun- | 

ctions in the Hartree-Pock procedure, obtained a value of -0.724. j 

The work of Neumann and Moskowitz"^^ has yielded a value of I 

—0.697. These authors employed an extended gaussian basis set. | 
Snyder and Basch^^ using a SOP wave function of double zeta i 

quality obtained a value of -0.3029. Comparing this value with | 

that of Nesbet brings out clearly the importance of inclusion of i 
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Table III. 10 

Atomic Orbitalwise Contribution to EPG in CO at '*'^0 Site 


Contribu- 
tions from 

CNI}0/2 

lEDO 

M-ITOO 

Oxygen : 




28 2s 

0.01345 

0.02287 ' 

-0.02269 

2p 

3.96680 

3.88786 

2.98350 

2p 

-4.20681 

-4.31376 

-3.62852 


-0,22656 

-0.40303 

-0,86791 

^^nuc * 

0.0 

0.0 

0.0 

Carbon ; 

2s 

-0.16080 

-0.15640 

-0.07963 

2p 

-0.00066 

-0.00300 

0.02554 

2p 

-0.33331 

-0.32192 

0.19097 


-0.49477 

-0.48132 

-0.24486 

^nuc » 

0.82550 

0.82550 

0.82550 

^total 

0. 10417 

-0.05885 

-0.28748 
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po iHTisation functions. Recently, Amos^^ made a SCF-C I calcu- 
lation oi the F'J* (j at the ^*^0 site in CO and studied its vihra— 
tiotuil do}u,-n<'i(‘ncu? also. The EFG value drastically changes 
from — o).6532 to -0.8436 with the inclusion of configuration in- 
teraction in the SOP MO wave function. The experimental value of 
—0.7168 due to lovas and Tiemann^^ ,is in good agreement with the 
rooont ah Initio results. At the semi-empirical level, ours is 
perhaps the first calculation that includes one- and two-centre 
intc/'iru. Is , The NN-INDO result is. very close to that of Snyder 
,'tnd Hanch iind is better than the CNDO/2 and INDO results. 

I 1 1 . 6 Sonic nomnii'-nta on the Role of Orbital Exponents in SPG 
Oalcu la t ions 

The tendency of NN-INDO to yield values of total ERG nume- 
rically higher than the other two NDO schemes is amply clear 
from the present work. This has its origin mainly in the use 
of AO’s with exponents derived from Burns rules rather than 
Slater rules. Burns rule AO's are generally more spread out 
than the STO’s. Whitehead^^ compares' the effect of Burns ex- 
ponents in a molecule such as IIT' with that of Slater exponents 
and the IIIjMO exponents, and finds that the Burns exponents tend 
to accumulate charge in the inter-nuclear region near R, giving 
rise to a numerically low value for the electronic field gra- 
d iont for the H site. Our results on linear molecules, also show 
that the NN-INDO method, using Burns exponents,does give too low 
an electronic field gradient at the H site/ however, the electronic 
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contribution at the other nuclear sites in the molecules 
exarain.-^d here is uniformly higher in M-IJTDO as compared to eiI)0/2 
and INDO, th(' solo exception being LiD. The explanation for 
this lies in the fact, that with the greater 'expansion' of the 
orbita Is ,thf? two-centre terms tend to become more important. 

111.7 Empirical Evaluations of EFG's Using Electronic Population s 
on the Atom with the Quadrupolar Nucleus 

Several empirical methods^® for the calculation of EFG's have 

their basin in the fact that the one-centre integrals are far 

larger than the two-centre or three-centre integrals. For atoms 

other than hydrogem, negl-'cting contribution from cores, and 

assuming that other atom contributions and the nuclear contribu- 

49 

tions ouncul each other, leads to ^ the equation . 

„ n„ + n,, , ^ 

q - 9° (n^ ^^“2 — (III. 5) 

where q” is the due to a single p^ electron at the site of 
the nucleus, and n^, n^, and n^ are the populations in the p^, Py 
and p, orbitals on the atom in question. This equation, or bhe 
related one obtained by multiplying throughout by eQ, forms bhe 
basis for the empirical evaluation of QCC ' s or EFG s. 

We shall briefly examine here the validity of this type of 

14 

approach, using our MO data on the molecules containing N 
and nuclei. Clearly, this approach is meaningless for 

deuterium couplingr. Two kinds of populations are generally 
employed In eqn. (III. 5). are the 'net' and 'gross' popu- 
lations^^. In connection with the Croo/2 and related 
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thero is thn question of the nature of the orbital basis, too. 

If tho GMDO function is interpreted as being over an OAO basis, 

then with Botsuyulcu^', „e have to use the elements of the charge 
denoit.v-bonrl orhor matrix over the deorthogonalised basis and 
U.S "till! popiild "fcions "th.'u.s ob"t8.iii6d in 1^11© o^bovs fomn'u.ls.v In 
principle, we have to take the gross population over the over- 
lapping basis for the calculations. More commonly, the net popu- 
laticn frrni the calculation over the orthogonal basis has been 
used. This is unsatisfactory, since net populations over OAO’ s 
have C('nip(U\.'ntn of other orbitals in them and the procedure 
effect tv < ly Introduces many-center contributions. q° in eqn. 

(in. 5) thus loses its atomic significance. 

Thero does not appear to be a unique value of q° for each 

52 o 

atom. Literature values of q for nitrogen, for example, | 

range from -7.4 ti.) -14.0. In the present calculations , we have ; 

adopted the q® values obtained f.rom the <(^l/r^^ values estimated 
by Barnes and Smith^^ from atomic spectral data. This -ckroice is 
made because it is independent of any reference compound. Our q ra- j 
lues havt) been computed using four different schemes over OAO basis 
functions using net and gross populations. The results are pre- i 
sontod in Table III. 11 for the three semi-empirical MO formalisms | 
choaon. On comparing the empirical q data with experiment , there | 
does n'.it appear to b‘=' a clear correlation. Indeed, recently I 

Barber at al^f , based on their ab initio SCF, calculations on 
some nitrogen containing molecules, have also noted that there xs 
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Table III. 11 

v'onu rlcfil eTaluation of BFG-'s in Some linear Molecules 

using Populations 

(nuoUus t)X interest is indicated by an asterisk) 


Pormali sm 

4 

Basi s 

Population 

type 

N=N* 

HC=N* 

1BG-C=N* 

0=0* 

CMDO/2 

I 

Net 

-0.620 

-0.726 

-0.624 

-0.394 



Gross 

-0.900 

-0.754 

-0.822 

-0 . 700 


1 1 

Net 

-0.297 

-0.243 

-0.147 

-0.036 



Gross 

-0.398 

-0.419 

-0.397 

-0.229 

IdDO 

X 

Net 

-0-541 

-0.648 

-0.560 

-0.498 



Gross 

-0.831 

-0.690 

-0.774 

-0.774 


1. 1 

Net 

-0.219 

-0.157 

-0.072 

-0.263 



Gross 

-0.287 

-0.332 

-0.328 

-0.337 

NH-IN.UO 

II 

Net 

-1 .283 

-1.341 

-1 . 207 

-1 .327 



Gross 

-1 .089 

-1 .007 

o 

o 

o 

T 

-1 .829 

Exptl. 



-1 .025 

-1 .131 

-0.942^°^ 

-0.728' 


§ I indicates without deorthogonalisation and II indicates with 
doorthufjonalisation 


(a) Sou Ku,r. 33 

(b) uoo lt.il*. 37 

(c) Sue itof, 38 

(d) Sou Kuf. 46 



no evidence for any correlation between EFG’ s and Mulliken 
populciion analysis indices. The situation might not be different 
in the cese of other nuclei as well. Hence, in our opinion, for 
a quantitative evaluation of the merits and demerits of the va- 
rious MO wave functions, with regard to the prediction of EPG-'s, 
empirical methods based on population are unsatisfactory. One 
must therefore, evaluate directly the expectation value of the 
WPG operator over the MO wave functions, as we have done here. 

In the next two chapters, we present our BPS results at the 
and ''”^0 sites in several non-linear molecules. 
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IV^JL gene ral 

Wlt'h the exception of syinmetric top molecules » the principal 
axis of the EPg tensor need not coincide with the molecular sy— 
mnietry axis of non-linear molecules. Therefore, it is necessary 
to evaluate the off-diagonal and diagonal elements of the EPg 
tensor. Now, by the Laplace equation two of the diagonal com- 
ponents uniquely determine the third, and as the EPg tensor is 
symmetric, only three off-diagonal elements are independent. 
Hence, a total of five EEG components are required to describe 
the tensor completely. The EPG tensor can then be diagonalised 
by a similarity transformation 

C q 0"“' = q^ (IV. 1 ) 

From the diagonal EPG tensor q^, the principal component and 
the asymmetry parameter are calculated. In order to define the 
orientation of the EPG principal axis system with respect to the 
molecule-fixed coordinate system we further need the direction 
cosines and these are obtained as elements of the transformation 
matrix C, 

IMiMlIlll 

0^^ = cos(Uj_, V^); (i,j=1,2,3). (IV.2) 

where u^ , Ug and u^ refer to the molecule -fixed coordinate 
system and v^ , Vg and v^ refer to the EPG principal axis 
system. It is customary to represent such a rotational trans- 
formation of coordinates by using Euler angles connecting the 
two systems. Even though, such a procedure is more elegant, in 
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this thesis vre have chosen to represent the EFG tensor components 
in terms of angles relative to the various bond directions; this 
method has the advantage of easier visualisation and inter- 
pretation from chemical point of view. 

We follow the convention of defining the EFG- principal com- 
ponents q^y and q^^ such that 



The molecules we have chosen for study in this chapter are 
the 1. acyclic molecules CD^, D^O, and h2GO and the 

cyclic molecules furan and pyrrole. All our EFG calculations 
reported here have been done with the experimental geometry, 
obtained from the literature, preferring microwave data where- 
ever available. The geometries used are indicated in Fig. IY.1 . 
The rosults of CKDO/ 2 , INDO and M-IIDO investigations on the 
BFG's on these molecules are presented and discussed in the 
following sections. 

IV.;J EFG *8 in Home Acyclic M olecules 

We shall now consider the field gradients at the D site 
in GD^, NDj, JDgO, N2D4 and D2CO, the N site in and ir2D^» and 
the 0 site in HgO and HgCO. 
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FlG»IV*l- GEOMETRY OF THE MOLECULES STUDIED 
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IV. 2* 1 Deuterium Field Gradients 

The.' rniiults obtained using the CNDO/2, IHDO and M-IEDO 
achci.ujo for the principal component and asymmetry parameter 
for tho B-sitos mentioned above are listed in Table IV. 1. As 
has been already observed in Chapter III, the values obtained 
for deuterium field gradients by the NN-INDO procedure are 
higher than that from the other two methods. 

Our Gl^DO/2 and INDO calculations agree with Bloor and 

1 2 
Maksie'e value of 0,593 and Barfield et al.*^ 0.409, respecti- 
vely. Arrlghini et al.^ using 39 STO's in the ab initio SCF MO 

A 

calculation obtained a value of 0.316, Later Pranchini et al. 

improved upon this calculation by using strongly orthogonal 

geminals and obtained a value of 0.324. Snyder and Basch^ (see 

also Snyde^r^) used STO's of double zeta quality in their . ■ 

ab initio calculations and obtained a value of 0.3362, Dixon, 

Claxton and Overill*^ optimised the C-H bond distance and the 

orbital exponent on D and reported values ranging from 0.3405 to 

0*304 for exponents ranging from 1.3 to 1.5 at an inter-nuclear 

8 

distance of 1*95 bohr. Bend azzoli, Dixon and Palmieri using 
18 STO's report a value of 0,549. Dixon^ using double zeta, and 
other gaussian wave functions reported values ranging from 0.306 
to 0,352. Eberhardt, Moccia and Zandomeneghi^^ using 39 STO's 
and configuration interaction involving all single and double 
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Table I?,! 




EPG‘s 

at D Sites 



Molecules 

CNI)0/2 

'^zz 

n 

INDO 

^zz 

T) 

HH-nroo 


0.39608 

0.0 

0.40869 

0.0 

0.53256 

0.0 

ND^ 

0.47630 

0.10 

0.46706 

0.10 

0.70421 

o.n 

DgO 

0,48897 

0,13 

0,46400 

0.14 

0.67370 

o.io 


0,45765 

0.19 

0.44086 

0.20 

0,61006 

0.24 

D 2 C 0 

0.36432 

0.009 

, 0.38208 

0.007 

0.49804 

0.004 
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excitod conf igurstions of "A” synaaetry obtained a value of 

0.5209* Thu molecular beam electric resonance ejqjieriments of 

.11 

Wofsy et al* »load to a deuterium wCC value of 192 KHz, 
corresponding to an EFG value of 0.286 a.u. The experimental 
asymmetry parameter value is zero and is in agreement with all 
the reported theoretical calculations. However, the experimental 
q[ value is lower than all the values theoretically calculated. 

The CWDO/2 and INDO values are lower than the M-INDO value. 

This foaturo of NW-INDO was also noticed earlier with other 
douturium BFG's (vide Chapter III). 

HD, 

.j 

Our CNDO/2 results agree substantially with Bloor and Haksic’s 
calculations . These authors report a value of 0.395 for 
and T) sf 0.10, Their value for the angle between the H-H bond 
direction and the principal axis of the EFG tensor is 2.9°. 

Our INDO calculation gives for this angle the value of 2.1°, 
and with the NN-INDO method the EFG principal axis almost 

2 

coincidoB with the N-H bond. The INDO results of Barfield et al. 

using the Slater orbital exponents, are q^^ = 0.464 and n = 0.12. 

Our INDO T) value is slightly lower (0.10). Barfield et al. 

obtain a n value of 0.08 using Hehre-Stewart-Pople orbital 

exponents, Franchini et al*^, using 29 STO functions in their 

ab initio calculations, obtained values of q = 0.4286 and 

h s* 0.152. Using 28 localised geminals these authors reported 

. 12 

values of q « 0.4254 and r) = 0. 149* Kari and Gsizmadia , using 
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a GTO basia»qucte values of q = 0.4331 and n = 0,133. With 
confi^qaraticm int..raction involving 918 configurations ,Eari 
and CslKmn.dia obtained q = 0.4363 and t) = 0.128. Snyder and 

K C 

Btisch (oec.‘ also Snyder ) obtained using double zeta basis 

values of q = 0.459 and t) = 0,127. Bendazzoli et al.^ quote 

values of q = 0,457 and r\ = 0.120, with their STO double zeta 

functions (See also Dixon^). Cook and Ebbing^ ^ with their 

'coro' mndol calculation obtained values of q = 0.4977 and 

■n =0,09. Goddard et al.^^, using basis sets which incorporated 

tho conjugate gradient optimised N-atom basis, obtained after 

optimisation and scaling a value not very different from that 

1 ? 

of Rari and Caizmadia . The experimental values due to 
Cederberg are q * 0,308 and n = 0.072. Reviewing all the 
calculations in the light of this data, we see that the q-value 
has boon overosti mated by semi -empirical as well as ab initio 
formalisms. As far as the M-INDO method is concerned it gives 
a higher value than CNDO/2 and INDO. 

D2O 

2 

Our rv.'sults on method and Barfield et al.’s 

value of q„„ = 0.474 and q = 0.14, agree fairly well. Bloor 
and Makeic^ do not quote any ChDO/2 values on D 2 O, tut using the 
solf-ccnsiotcmt charge method they report a Igz ‘^alue of 0.534 
and h value of 0,10. The principal axis obtained by our 
CNDO/2 calculation lies in the molecular plane and makes an 
angle of 2,75® with the O-P bond direction. In the INDO 
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calcu Inti ono the value of this angle is 2,4°. in the IM-IKDO 

tho' nxia practically coincides with the 0-H bond. Kern 
1 6 

and W.Mtohr, Mtudi(,d the effect of vibrational corrections to 
tho q in wati-.i . iranchini, Moccia and Zandoneneghi^ report 
using 27 i.}TO * s q = 0,5413» h = 0,1183; using strongly ortho- 
gonal geminnl basis they roported values 0.5520 and 0,1145 
rc-opectivoly for q and p . Brmler and Kern"''^ with a GIO 
baoiB t)btai!\i;d q ss 0,502 and p =0,138 with vibrational corre- 
ctions nml q = 0.564 and p = 0.127 without vibrational corre- 

1 fi 

ctione, Krohn, Ermlor end Kern report q = O .569 and p= 0.129 

8 

with vibrntionnl corrections, Bendazzoli, Dizon and Palmieri 

using STO double zwta basis report q = 0.576 and p= 0,111; with 

OTO basis their values are 0.547 and 0.115 respectively (See 

also OiKOn^). Schaefer has discussed some of the ab initio 

calculntlims on EFG's at the deuterium site in water in terms of 

5 

the quality of the wave functions. Snyder and Basch report EEG 
components , from their double zeta quality STO MO wave functions, 
which load to q^^ and p values, 0,561 6 and 0.115 respectively. 

pA 

Rosenborg and Shavitt obtain q = 0 . 5625 and p = 0,115 using a 
doubii? Zijta bt^ais with polarisation. On the experimental side, 
Bluy , Verhoeven and Bymanus^'^ have reported values of 0,458 
and 0,14 for q and p respectively. It is gratifying to see 
that b'-ith GNiJO/2 and IKDO values are in good agreement with the 
cxporiin<intal values. The M-INBO again gives too large a q 
valuQ but a reasonable value for p . 
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^2 4 

Our roaults are based on a conformation where the two HD2 

groups arc staggered with respect to each other. Snyder and 
5 6 

Basch (bee also bnyder ) used a different conformation in' 
their SCP MO calculations and obtained q values of 0.4274 and 0.42 
for the two deuterium eltee and t) values of 0.104 and 0.141 res- 

pectivuly. Comparing oui' q values with those of Snyder and 
Baach we see that the CNDO/2 and INDO values are in reasonable 
ngroomcnt while the BN-IWDO gives too high a value. With our 
geometry wo find that the diagonalised BPG tensor has a prin- 
cipal axis, which although does not coincide exactly with the 
N-D bond lias within 1° of it. To our knowledge, no experimental 
data on tho KFG'e at the deuterium sites are available. 

BgCO 

CXir CN]}0/2 values on DgCO agree well with Bloor and Maksic’s 
CI^BO/2 value of O.364 for q and 0.008 for n . Barfield et al.’s 
IWBO rfjsults^ on this molecule are based on the Hehre-Stewart- 
i’ople orbital exponents. The values obtained by Barfield et al. 
in this manner are q = 0.575 and p = 0.01, The agreement with 
our IBBO values using Slater exponents is reasonable. .’'ITeumann 
and Mnekowitz^^ using contracted geminal functions obtain 
q r= 0,262 and n = 0,157, Snyder and Basch^ using SOP MO 
wave functions of double zet a quality have obtained q = 0,286 
and T) = 0,007, Garrison et al.^^ have reported components of q 
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Table IV*2 
EPG’s at N Sites 


Moleculoo 

CNDO/2 


iroo 


m-imx) 




n 

^zz 

T) 

q T) 


-1.45798 

0.0 

“1.59082 

0.0 

-1.460140 0.0 


0.95222 

0,026 

0.88208 

0.163 

-1.17287 0.214 
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for a ^/nxjGSjan basis. On diagonal! sati on, the values we obtain 
ero q “ 0*21964 and q = 0.184. . This r\ value appears to be 
too large. Experimentally Plygare^^ obtained a value of 0.2622 
for the Hold gradient component along the G-H bond axis, with 
T] ■- 0.02 whereas Tucker and Tomasevich^^ obtain values of 
q = 0.2426 and h = 0.04. 

In our calculations the principal EPG axis has been deter- 
mintui in nil the three schemes. In the CNDO/2 scheme we get for 
thu nnglf) butwcon the C-H bond and the BEG principal axis a 
vnluo of 1,3“. In IJMDO, this angle is changed to 1.6°, while 
JNN-INDO puts tho angle as 3*6°, in the opposite sense. It is 
worth noting that all the three semi -empirical MO formalisms 
used hero yield higher q values and very low n values in com- 
parison with experiment, 

IV. 2. 2 Nitrogen Field Gradien ts 

BFG's at the W site in and have been studied in 

the present work. Our results for these are presented in 
Table IV, 2. 

NHx 

Thu calculated principal BFG axis for the N-centre in 

ammonia from GWDO/2, I EDO and M-IKDO formalisms is found to 

26 

coincldo with tho molecular symmetry axis. O'Konski and Ha 
invootigatud the BFG's at the N site in using nine and 
twelve gaussian lobe group orbitals and obtained theoretic al.q 
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C)i •*‘1 *959 Slid *•1*045^ ^^sspoctively* Eletr et 
iiiv^ 'U ti, /V 1 1 i:d lhc‘ c’l i uctiveneas of one-centre expansion method 
and ohb«iint‘d ruoults in satisfactory agreement with multi- 
Ci.'nii*.* w«iV‘t‘ iunctions* Ditchfield et al*^^ report values 
rf-in^ln^*, Iroin -1 *277 to -1*732 for ST0-n& functions with n ranging 
from 2 to 6. Olympia^^ examined Configuration Interaction (Cl) 
otiu'CtiT on th.j iil'G at the centre in using fixed molecular 
gt()mi.’ti\y nnd orbital parameters with a series of one-centre expan- 
oioruu With a single configuration the q was found to be -0.8133, 
with 01 the value became -0.92951 . Dewar et al.^^ calculated the 
KlU in Wii.. but those authors approximated the two- and three- 
coritre Integrals and hence their results are not strictly com- 

9 

parabli.) with oither our results or ab initio results. Dixon 
calculatod the EFU in NH^, using BTO double zeta wave functions as 
well as vf.^rlcus gauasian expansions for STO basis, and obtained 
a value of -1,315 • for the double zeta function and values 

ranging from -0,602 to -1,281 for the various gaussian expansions. 
D.A. Zhogolov et al.^^ used a mixed basis of STO’s and gaussians 
con trod on E atoms in their ab initio calculations. Their value 
for this basis was -1,76, agreeing with the pure STO result of 
-1.7/1. harbor, Hayne and Hinchliffe^^ have recently obtained 
an ab initio value of -1,0471* Experimentally the value obtained 
by Kukolich'^'^ » by molecular beam method for q>is -0.877. The 
ab initio ri'sults cover a wide range of q values while the semi- 
empirical results are higher than the experimental value. The 
NN-IEDO result is not much different In this case from the CNDO/2 
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find 1. Mix) rt.,'Gult 3 , This situation is to be contrasted with 
thf; diiuturium case, 

3'4 

Arri;':hinl et al, ^ using gaussian lobe orbitals calculated 

a. q value of -2.390 and a t] value of 0.561. Kochanski, Lehn 

and Lovy-^^ obtained a q value of -1.56 for the nitrogen field 

gradient in along with a n value of 0.3846. Snyder and 

Sfisch's^ calculation with their double zeta wave functions lead 

to a q value of -1 .2793 and a q value of 0.5351. The recent 

32 

rctjults of Barber;. Hayne and HinOhliffe are q =-l.29 and 

36 

n c 0,778. Harmony and Baron from their microwave studies 
roportfid a value of q =-0.9328 and a n value of 0.064. Coming 

now to our computed q values, ONDO and IHDO values have the 

positive sign whereas M-INDO value has the negative sign, and 
is a definite,' improvement. Our p values are, however, closer 
to tho microwave data. 

IV. 2. 3 O xyg en Field Gradients . 

Wo have examined the ElG's at the ^'^O sites in water and 
formaldehyde. The results are presented in Table IV, 3- 

^ I 

The principal component of the EFG- tensor has a positive I 

sign in the Nj) 1-IMH0 method, whereas it has a negative sign in j 

CHDO/2 and IHDO. Also, the direction of the principal component | 
is found to be perpendicular to the molecular plane in the OT-INBq 
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Table IV. 3 
EPG' s at 0 Sites 


Moleoulo. 


IIDO 


M-IUDO 

^ZZ 

IQ 


T) 

^zz ^ 

H^O -2.93176 

0.62 

-2.92281 

0.67 

2.26314 0.37 

HgCO -2.64430 

0.07 

-2.50956 

0.01 

-1.03152 0.85 
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and parallul to the H-H direction in CNDO/2 and lUDO. Snyder 

and Bnsch’o roported theoretical EPG component values^ lead 
1 7 

to thu C) princijal axis being perpendicular to the molecular 

plnno, but their q value is negative (- 2.0487). Aung, Pitzer 
37 

and Chan obtained, using a STO minimal basis set, a q value 

of -2»6l6. Korn and Matcha obtained, using a GIO basis, the 

values of q « -1 .922 and h = 0,796 with vibrational corrections. 
3B 

IXinning using a GTO basis set got q values of -1,9711 and 

-1 .9245 uuing different basis set contractions, Gornostansky 

39 

and Korn , using a BTO minimal basis set obtained a q value of 

-2,62 and q value of 0,78. With STO double zeta basis set 

Bendazis'.oli , Dixon and Palmierl^get a q value of -2.288 along 

with a q value of 0,7928, Some of the ab initio results have 

been aummarisud by Schaefer ^ in his book. Por some recent 

9 

ab initio results we may refer the reader to Dixon . A more 
recent ab initio calculation on water is that of Rosenberg 
and Shavitt^^, Their values are q = -1,729 and q = 0.799 . 
Using rriasor beam spectroscopy Verhoeven^^et al. obtained an 
experimental value of -1.6463 for q and a value of 0.75 for q . 

At the a omi -empirical level ours appears , to be the first 
calculation with the inclusion of all integrals. The CNDO/2 and 
IhDO results predict correctly the sign of q. They also yield 
a fairly reasonable estimate of q . The performance of M-IFDO 
however, is rather disappointing. Whether this is a consequence 
of the paranetrisation or is due to a deep-rooted cause is not 

clear at this stage. 
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li^CO 

In all our three calculations based on GNDO/2, INDO 
and NN-INUO the HIPG principal component direction is found 
to be in the molecular plane and perpendicular to the C = 0 
bond , 

00 

Neumann and Moskowitz using contracted gaussian basis 
sots obtained a value of -2.2705 for q and 0.644 for p, 
with the principal axis perpendicular to the C = 0 bond. They 
alao commented that the earlier work of Poster and Boys^^ ■. 
was In error mainly due to the error in geometry. We have also 
noted that our calculated results are quite sensitive to 

5 

ohangos in assumed geometry. Snyder and Basch obtained a q 
value of -2*6124 and a ri value of 0.467» with the axis of 
the principal component of BPG lying in the molecular plane 
and perpendicular to the G = 0 bond. Ditchfield, Miller and 
Poplu^®, using various STO-nG and STO bases obtained values for 
q ranging from -2*194 to -2.9294, the value for the ST0-6G cal- 
culation agreeing with the STO result. Their p values range 
from 0.640 to 0.837, the last referring to the STO calculation. 
Garrison at al.^^ report a value of -2,302' for q and 0.624 for 
p • Flygare and Lowe^^ have obtained from their high resolu- 
tion microwave, spectroscopic studies a q value of -2.001 and 
a p value of 0.695. On reviewing the semi -empirical results 
we see that the q values of ONDO/2 and INDO appear to he better 
than the NN-INBO value. However, the t) value from M-INDO is 

superior. 
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J,n_ g ome Oyclic Mo le cules 

■do huvu attlocted for our study two typical hetero-cycles 
ono cnntn i nl n/'^ nitrogen and the other containing oxygen, viz., 
pyrrolu and furan. The EPO results on these molecules are 
summarised in Tables IV. 4 and IV. 5 respectively. 

I'.vrrole 

BKU'b at tho nitrogen site in pyrrole have been studied 
by Kochannkl , Lehn and Levy using contracted gaussian fun- 
ctions. Their q value for the nitrogen is -1.40 which 
corrifUi^ondn to n QCC of -6*348 MHz using the conversion factor 
given in Appendix II. The QCC value given by Kochanski et al. 
is -*5.24 Miz since they use a different conversion factor based 
on the older Q value of Their calculated n value is 

0.109. Von Niossen et al.'^^ report, from theif ab initio studyr 
a q valuo of -0.8491 and a q of 0.07. The experimental QCC 
value from microwave studies^^ is -2.66 MHz which corresponds 
to a q value of -0.587. The q value reported in this study 
is 0,09. The performance of all the three semi-empirical 
formalisms is poor. This is not surprising in view of the 
fact thot oven with the ab initio approach von Miessen et al. 
woru not able to reproduce the q value. For the hydrogen 
attached to the N'r von Hi ess en et al. obtained the values 
q = 0.3691 , and h • 0.155; the EFG principal axis for this 
hydrogen lies along the N-H bond. The a-hydrogen has q-0.3625 
and Ti 0.081 and the principal EFG component makes an angle 
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Table I¥.4 

EPG Parameters for Pyrrole 


Site 

CNl;0/2 

q 

T1 

INDO ■ 

^zz 

r] 


) 

r\ 

N 

-1.58129 

0.409 

-1.57936 

0.399 

0.59388 

0.968 

H(N) 

0,37913 

0.122 

0,38686 

0.120 

0.52791 

0.056 


0,41413 

0.033 

0.42389 

0.039 

0.45899 

0.029 

Ha 

0.40353 

0.018 

0.41204 

0.024 

0.44130 

0.004 
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Table 17.5. 

BPG Parameters for Puran 


Site 

nNl)0/2 

Hi <f 

r\ 

IITDO 

^zz 

T1 

M»IND0 

^zz 

rt 

0 

1.62241 

0,979 

1,71467 

0.750 

-0.92936 

0.990 

a 

0.43470 

0.024 

0.44731 

0.030 

0.48349 

0.017 

U 

0.42539 

0.005 

0.43529 

0.019 

0.45743 

0.013 
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of 0,46° with the Cq.-H bond. For the p-hydrogen von Kiessen 
et al. calculate q = 0,3656, rj = 0.058 and the angle between 
SPG principal component and Cp-H bond as 0.12°, In our cal- 
culations the principal component of the EPG tensor at the 
hydrogen attached to the a carbon makes angles of 1.6°, 3*5° 
and -4-5° respectively, with the Ga-H bond in CMDO/2, IlDO and 
NM-IMDO formalisms, while that of the hydrogen attached to the 
P-carb on makes angles of 5°, 5° and 5.3° respectively with the 
Gp-H bond in the above three formalisms. As far as the 

orientation of the principal EI'G component at the K site is 
concerned, it lies perpendicular to the molecular plane in all 
our calculations. 

The agreement between the results of the semi -empirical 

formalism and von Messen et al.’s ab initio a.pproach in the 

case of deuterium EPG's is good and is surprising, in view of 

1 4 

the failure of the former methods in the case of N. 

Fu r an 

Von Niessen et al. have been quoted by Moccia and 
45 

Zandomeneghi as having obtained from their ab initio calcula- 
tions the values of 0.3584 and 0.3619 for the EFG at and Dp, 
with asymmetry parameters of 0.073 and 0.053 respectively. For 
the oxygen centre, to our knowledge, there is only one earlier 

AC 

report*’ of a calculation but the EFG values from this work not 
readily available to us for comparisons. The deuterium field' 



gradients from semi -empirical methods appear to he consistent » 
viewed in the light of the ah initio results of Yon Niessen et al. 
However, the 0 results do not show the same consistency. Ihe 
sign of q obtained in HN-IImDO for the oxygen site is different 
from that of CHDO/2 and INDO. Our .calculated EPG tensors for 
deuterium sites in this molecule have their principal component 
practically in the same direction as that of the corresponding 
C' D hond. 

Before concluding this chapter it may he remarked that while 
the deuterium EI’G's .appear to he amenable to theoretical treatment 
the nitrogen and oxygen EEG's offer a challenge of a high order. 

We shall now examine in the next chapter a few mere systems 
containing deuterium, nitrogen, oxygen as well as boron. 
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CHAPTER 7 

EFG's AT E,B,N AND 0 SITES IN SOME MOIBCUIES 



V. 1 General 


In Chapter III we saw that the EPS at a given nuclear 
site is . largely determined by contributions from that atom 
and its near neighbours. This qualitative generalisation 
fonae the basis for a search for chemically significant 
regularities among the EPS' s in a number of related systems. 

With this background in mind we shall examine, in this chapter, 
the EPG’s at the various quadrupolar nuclear sites in the 
following molecules. These are : CD^CN, CD^ND2, KJOND2, 

D202,D01I0, B2Dg,B^lT^Dg, BF^, N25'2» ^'30 and P2O2. The geometries 
of the molecules studied are indicated in Figure V, 1 , In cases 
where the EFG's are non-axial we report the asymmetry parameter 
values as well, 

V .2 EFG»s at D-Sites 

We shall first consider the deuterium BPS* s in the series 
CD^,G2D^, C2B2 (see Table 7 . 1 ). The essential distinction as 
far as the C-D moiety is concerned is the state of hybridisation 
of the carbon. Methane and acetylene have already been discu- 
ssed in the previous chapters, but the values are reproduced 
here to enable an easy comparison along with our ethylene 
results in Table 7 . 1 . Bloor and Maksic^ obtained for ethylene, 
by GNDO /2 method, the values q = 0.413 and n = 0 . 029 , which are 
somewhat different from ours,* the difference is probably due to 
the difference in geometries used. Barfield et al. ,using the 
IlTDO formalism with Hehre-Stewart-Pople orbital exponents 
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FIG. V. I (a) . GEOMETRY OF THE MOLECULES STUOtEO 



118 



FIG. V.l (b). GEOMETRY OF THE MOtECULES STUDIED 
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Table IT.l 

EFG's at D Sites in and 


Molecule 

Method 

<1 

n 


croo/2 

Q. 39608 

0.0 

CD4 

INDO 

0.40869 

0.0 


NN-INDO 

0.53256 

0.0 


CNDO/2 

0.46846 

0.027 

°2®4 

lUDO 

0.47791 . 

0.033 


ra-INDO 

0.59901 

0.027 


OEDO/2 

0.63291 

0.0 

CjDj 

lEDO 

0.62369 

0.0 


NN-INDO 

0.53022 

0.0 


t20 


obtained a q value of 0.4468 with ti = 0,04 for ethylene. Gal- 
culations at the ab initio level for this molecule have been 
done by Harrison^ and Schulman et al.^» yielding values of 
q = 0.3780 and q = 0.3188 respectively. The latter authors 
obtained a t) value of 0.05. Snyder^ has reported a value of 
q = 0.344 ^nd rj = 0.039 using double zeta basis set for 
ethylene. 

It is seen clearly that CNDO/2 and IFDO methods give q 
values that decrease along the series €2^2 ^ ^2^4 ^ ^^4* 

5 

same trend can be seen from the ab initio results of Snyder . 

On the other hand, the NN-INDO results do not show this trend. 

The exponents used on the carbon atom in the !S[lT--IErD0 method are 

13 

the Burns exponents which correspond to a valence state s p . 

It seems likely that due to this choice of exponents, the deu- 
terium q values do not reflect the above trend in the HH’-INDO 
calculations. The h values for ethylene from UH-IKIX) agree': 
fairly well with that from other methods. We may note here that 
Rinne, Depireux and Duchesne^ had pointed out earlier the depen- 
dence of the EPG- at the deuterium site in the C-D bond on the 
state of hybridisation of the carbon. The trend ^2^4^ ^^4 

noticed here with CNDO/2 and IRDO formalisms is in conformity 
with the work of these authors. 

Table '7.2 presents our results on the molecules CB^CR, 
CB^IIH 2 and DCOIS 2 . Bloor and Maksic^ report a value of 
q = 0.378 and h = 0 by the CNDO/2 method for CB^N. Onr ClIBG/2 
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Table T.2 


EPG’s at D Sites in GD^CF, CD^KH2»^0HH2 


Molecule 

Method 

q 

D 


CNDO 

0.35576 

0.003 

CD^CN 

IUDO 

0.36578 

0.000 


NU-INBO 

0.48090 

0.020 


CNDO 

0.39520 

0.024 

CD^EE2 

INDO 

0.41039 

0.022 


M-INIX) 

0.52945 

0.008 


CWDO 

0.43462 

0.016 

]DCOm2 

INDO 

0.45222 

0.009 


NN-INDO 

0.52734 

0.010 


t22 


values are q = 0.356 and t) = 0.003. The non-zero q in our 

P 

calculations is to he noted. Barfield et al. were able to 

carry out their calculations on GB,C1T with the inclusion of 

3 

one- and two-centre integrals only. Their values at this level 

of approximation with Slater and Hehre-Stewart-Pople e 2 i>onents 

are 0.499 ^nd 0,437 respectively with p = 0.05. Our IlsIBO values 

with the inclusion of all integrals and with Slater exponents 

are q = 0.366 and p = 0. The difference between q = 0.499 

and q =0.366 viz. ,0.133 is clearly seen to be the effect of 

the three-centre integrals. The NN-INBO method yields a higher 

q value in comparison with CirDO/2 and INDO as in the other 

deuterium calculations. The asymmetry parameter value is also 
15 

higher. Snyder has reported values of q = 0.309 and “n = 0,049 

7 

from his ab initio studies. Eukolich et al. ' have experimentally 
determined the q value to be 0.249. All the theoretical q values 
are seen to be higher. 

For CD^1IE[2 Harrison obtains, using gaussian lobe functions, 
ab initio q values ranging from 0.3552 to 0.3596 depending on the 
conformation chosen. Our CIH)0/2 and INDO values are in good 
agreement with these ab initio values. The HN-IKDO values are 
again on the higher side. 

For DCOHHg , Snyder^ obtains q = 0.307 and q = 0.03 from his 
ab initio calculations. Our semi-empirical q values are all 
higher while p values are lower. 
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Comparing CD^CN and CD^1IH2 , we see that the q Talue for 
the D-site in CD^N is lower than that in CD^NH 2 . This trend • 
is maintained in all the three semi-empirical formalisms. 

We may attribute this effect to the electron^ativity of the 
substituent on the carbon to which the deuterium is attached. 

It may be recalled that in Chapter III we had analysed the 
EP6 at D-site in DON and DGCD atomwise and found that the 
major contribution, to the EFG at the deuterium site comes 
from the carbon attached to the deuterium of interest. Since 
this is a two-centre contribution it depends on the electron 
density on that carbon atom. With a reduction in the density 
at 'this carbon it is reasonable to ezpect a reduction in the 
electronic contribution to the SPG with a consequent increase 
in the total EPG. This conclusion is further supported by the 
trend in the q values of IX 3 ONH 2 and D2^0 (0.4546 and 0.5645 
respectively, from CNBO/2 ealculations, for example). 

Results for the BOOB and BONO molecules are presented in 
Table V,3. Snyder^ reports the values q = 0.574 and n = 0.146 
for BOOB, The values obtained by us from the semi-empirical 
formalisms are in fair agreement with regard to both q and ri. 

Our calculations on BONO are, we believe, the first to be re- 
ported* We find that the SPG at B is greater in BOOB than in 
BONO consistently in all the three formalisms. The electro- 
negativity effect is flearly reflected here too. 

Results for the B-sites in CH 5 NB 2 , HCONB 2 and the B attached 
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Table V.3 

EFGr’ s at D Sites in ^ 2^2 


Molecule 


D2O2 



DONO 


Method ■ q. 

CMDO/2 0.50794 0.130 

IMDO 0.50942 • 0.267 

NN-Iimo 0,69867 0.109 

CNDG/2 0.43034 0.112 

BIDO 0.40728 0.120 


0.56637 


NIT-IITDO 


0.086 
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to ¥ in borazine are presented in Table ¥.4. Between CH^KD 2 
and HC 0 ¥D 2 , the differences in the q values are small, being 
in a range that could be accounted for by simple conformational 
differences alone. The effect of the farther atom is clearly 
seen to be marginal. Harrison has obtained deuterium q values 
ranging from 0.4598 to 0.4607 for CH^HHD (depending on the 
conformation) from his ab initio calculations using gaussian 
lobe functions. Our values from 010)0/2 and IHDO are in 
excellent agreement with these values. To our knowledge no 
other EFG calculations on borazine are available. The q and q 
values obtained by us for the deuterium attached to the nitro- 
gen in this molecule are similar to those obtained earlier for 
the D-N moiety in other molecules. • 

Our results on the EPG at B-sites in diborane and borazine 
(D bonded to B) are presented in Table ¥.5. A remarkable fea- 
ture of the EFG’s at the deuterium sites in diborane is the 
fact that the terminal D-site has a positive q value while the 
bridge D-site has a negative q value. This feature was first 

5 

noticed by Snyder , who made ab initio calculations on this 

2 

molecule. Later, Barfield et al. ^ who performed I¥D0 calcula- 
tions with the inclusion of all Integra Is, also noticed this and 
have commented on it. Our CHDO/2, INDO and M-IHDO results all , 
confirm this feature. The asymmetry parameter predicted for 
the bridge D-site is higher than that for the terminal D-site. . 

¥0 experimental confirmation of these calculations are available 
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Table V.4 

EPff's at D Sites in CH^HL 2 ,HC 0 ND 2 and 

Borazine 


Molecule 

Method 

q 

h 


C 10)0/2 

0.46674 

0.096 

CH^ND 2 

INDO 

0.45722 

0.106 


lOr-IlTDO 

0.67585 

0.106 


CBDO/2 

0.45529 

0.150 

HCOIID 2 

INDO 

0.45961 

0.148 


M-IFDO 

0.67698 

0.118 


Cm) 0/2 

0.46204 

0.042 

Borazine 

INDO 

0.46757 

0.044 

(B 3 H 5 N 3 D 5 ) 

NN-INDO 

0.59184 

0.025 
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Table Y. 5 

EPG's at D Sites in Diborane and Borazine 


Molecule/ 

ITucleus 

Method 

a 

h 


CNDO/2 

0» 30627 

0.033 

^2^6/^term 

IHDO 

0.31430 

0.027 


M-IEDO 

0.32843 

0.003 


CITDO/2 

-0.36880 

0.252 

^2^6^ ^bridge 

ITOO 

-0.37186 

0.255 


NN-IFDO 

-0.45930 

0.193 


CNDO/2 

0.30199 

0.015 

Borazine/ D 

lEDO 

0,30797 

0.007 


M-INDO 

0.20590 

0.016 
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at present. The q values from semi-empirical wave functions 
are all higher than those from the ab initio calculations of 
Snyder. The q values for the D-site obtained in borazine and 
the terminal deuterium in 82 !)^ from the CNBO/2 and IKDO for- 
malisms are close to each other, reflecting the sp -hybrid 
character of the B orbital bonded to D in both cases. However, 
FH-INDO gives significantly lower q values for borazine. Con- 
sidering the fact that the NF-INDO method yielded uniformly 
larger q values at the D-sites in other cases, the situation in 
borazine appears anomalous. 

V.3 EFS’s at B-Sites 

Our results on the B centres in B 2 Hg,borazine and are 
presented in Table V. 6 . The q value for the B-site in diborane 
from the ab initio calculations of Snyder and Basch is -0.3074. 
The T) value obtained by them is 0.844. All our semi-empirical 
calculations yield q values which are lower in magnitude although 
their sign is in agreement. The asymmetry parameter values 
obtained from our calculations are also lower. Since no other 
calculations on the BFG-'s at the B site in borazine are available 
for comparison we shall not discuss them further here. In the 

Q 

case of Snyder and Basch report a q value of 0.2218 and an 

asymmetry parameter value zero for the B site. Our semi- 
empirical calculations using the Cl!nX)/2 and lEDO formalisms 
yield q values which are an order of magnitude larger. The 
NN-IHDO results are not available at the time of writing of this 
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Table 7.6 


EPG' s 

at B Sites in Biborane, Borazine and 

Boron Trifluoride 


Molecule 

Method 

q 

ri 


CBIX)/2 

-0.17509 

0.433 


HTDO 

-0.17037 

0.311 


M-INDO 

-0.19711 . 

0.064 


CNDO/2 

0.11211 

0.332 

Borazine 

IIDO 

0.08346 

0.562 


ra-iNix) 

-0.16135 

0.713 


CNI)0/2 

3.39730 

0.003 

BP 3 

INDO 

3.41369 

0.002 


NR-INDO 


— 
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thesis. However, judging from the performance of the semi- 
empirical MO wave functions , against the back drop of the ab 
initio wave functions ,for the prediction of EPO’s at boron 
sites we feel that the parametrisation for the boron atom in 
these formalisms would require further optimisation. 

V.4 EFG»s at H-Sites 

Our results obtained on the systems and 

CH^CEf are tabulated in Table 7.7. The conformation chosen 
for CH^KH2 staggered conformation, ClIDO/2 and IHBO 

results show an increase of the q (numerical) value from 
CH^NHg to NC0KH2» but NN-IHDO shows the reverse trend. 

For initio results of O’Ebnski and Ha^ 

as well as those of Arrighini et al,^*^ agree with each other. 

Their q value is -1,918 and the asymmetry parameter value is 

0.177. Our semi-empirical results are in reasonable agreement 

with these. Our results for formamide also appear to be in 

8 

reasonable agreement with those of Snyder and Basch . For 

CH^CN,our 01100/2 and lEDO results yield low positive q values, 

while the NH-IHDO yields a negative q value. Ab initio results 

ft 

of Snyder and Basch , as well as the recent ones of Barber 

et al. ^ ^ , point to a q value of -0,704 and -1.009 respectively. 

12 

The experimental value quoted in the paper of White and Drago 
is -0.9283* The M-IHDO seems to fare better than CfEDO/B and 
INDO in this case. 



13t 


Table V..7 . ■ 

SFG’s at K Sites in HCONH 2 and CH^N 


Molecule 

Method 

q 

9 


CNDO/2 

-1.42210 

0.288 

CH^KH2 

INDO 

-1*55330 

0.265 


ra-INDO 

-1.32682 

0.367 


ClTDO/2 

-1.91015 

0.160 

HCOIIH 2 

INDO 

-1.89174 

0.T58 


OT-INDO 

-1.23845 

0.161 


croo/2 

O.Q1951 . 

0.000 

CH^CN 

mo 

0.10415 

0.000 


M-INDO 

-0.39622 

0.000 
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Our results on HONO, 1T2J'2 a-ud ‘borazine are tabulated in- 
Table V.8. The field gradients in HONO by the three methods 
are comparable. . A similar situation is obtained for ^2^2* 

O 

The results of Snyder and Basch lead to a q value of -1.012 
and ri value of 0.47 in 1T25'2‘ '^^^6 initio q value is in good 

agreement with the IfN-IlirDO result* In borazine, the EFG’ s at 
N-site obtained from both GinX)/2 and INDO methods are compara- 
ble but the NN-INDO q value is nearly halved. The t) value 
obtained f rom , NN-IUDO is also smaller. 

V.5 EFG'^s at 0-Sites 

The oxygen field gradient parameters for the molecules 
^ 2 ^ 2 * ^ 2 ^ ^2^2 presented in Table V, 9. All the EPS 

calculations point to uniformly large T) values in these 
molecules. The ITIT-IFDO method gives positive q values in all 
these three cases whereas the Cin)0/2 and INDO methods yield 
negative q values. Whether this result is due to the differen- 
tiation of the p , p and p orbitals in the NN-INDO method or 
•^x' 

is due to Burns exponents is not clear, further work is 
required in order to understand the situation better. 

The q values for oxygen in the molecules HCOKH 2 HONO 

(vide Table V.IO) obtained from ClJDO/2 and INDO methods are ^ 
comparable. The high asymmetry parameter values for oxygen in 
HCONH 2 and the oxygen bonded to the hydrogen in HONO from CNDO/2 
and INDO formalisms is to be noted. On the other hand, the 
asymmetry parameter values for the terminal oxygen in HONO are 
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Table V,8 

EFG’s at F Sites in HONO, ^2^2 Borazine 


Molecule 

Method 

q 

T1 


CNDO/2 

-1.01908 

0.172 

HONO 

INBO 

-0.94938 

0.257 


NN-INDO 

-1.10890 

0.252 


CNBO/2 

-0.81625 

0.905 

^2^2 

iimo 

-0,80920 

0.886 


NN-INBO 

-1.15210 

0,803 


GBBO/2 

0.91661 

0.S06 

Borazine 

mo 

0.91360 

0.838 


M-INBO 

0.44907 

0.159 
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Table V.9 

BFG's at 0 Site in H 2 O 2 , ^ 2 ^ ^2^2 


Molecule 

Method 

1 

T1 


CIffDO/2 

-3.14657 

0 . 319 . 

H 2 O 2 

INDO 

-3.17426 

0.347 


NH-INDO 

2.13108 

0.588 


croo/2 

-3.58430 

0.338 

F 2 O 

INDO 

-3,69543 

0.297 


NN-IMDO 

2.56464 

0.795 


CIDO/2 

-3.12794 

0.438 

F 2 O 2 

INDO 

-3.07326 

0.419 


NN-INDO 

2.46060 

0.521 
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Tatle ¥.10 

EPG’s at 0 Sites in HCOKHg and HONO 

Molecule / Method q t] 

Nucleus 


HGONH2/O 


HONO/0* 


H0N0*/0<^ 


CNDO/2 -2.46681 0.741 
INDO -2.33154 0.786 
ra-INDO -0.96476 0.197 


CNDO/2 -2.73100 0.593 
INDO -2.62114 0.709 
NN-INDO 1.85059 0.180 


011)0/2 -2.97278 0.178 

INDO -2.85957 0.272 

-1.49687 0.250 


m-INDO 
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lower. The NN-IITDO procedure yields q values of opposite 
signs for the two oxygen sites in HOFO. As noted earlier, 

CITDG/2 and IBDO procedures yield results which are close to 
each other for oxygen EPG-'s, while the NIT— IlTDO procedure yields 
results which differ markedly. This implies that the electronic 
contribution to the EPG is estimated to be very much different 
in the NN-INDO procedure as compared to the other two procedures* 
A better picture can be obtained only after detailed ab initio 
as well as gas-phase experimental EPG results in this class 
of molecules become available. 
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COKCinDIITG REM^S 

We have examined in this thesis ,the BPG’ s at the D, Li, 

B,1T and 0 sites in a variety of molecules using the three 
semi-empirical MO formalisms: CinX)/2»INDO and OT-INDO. For 
the deuterium case, since the major contribution to the EFG 
comes from the atom to which it is bonded, the orbital 
exponents for the neighbour atom play a very important role. 
Since the use of Burns exponents results in a more diffuse 
orbital, the electronic contribution to the EFG at the D— site 
decreases, and consequently the total EFG increases in com- 
parison with the EFG values obtained with the use of slatei 
exponents. Using experimental EPG values, at the deuterium 
sites, it appears that it would be possible to optimise the 
orbital exponents in semi-empirical formalisms. The EFG at 
the Li site in LiL calculated by the NN-IKDO method is also 
encouraging. On the other hand, the EFG* s at B-and 0 -sites from 
the M-INDO appear to be anomalous. This clearly indicates the 
need for reparametrisation of the EN-IELO method. The hope 
that a formalism, .as the M-IEDO, which differentiates bet- 
ween, p , p and p orbitals would succeed in predicting the 
X y z 

asymmetry parameter values correctly has been belied. The^^^^ ^ ^ 
reason for this failure appears to lie in the fact that the 
behaviour of the p-orbitals near the nucleus is yet to be 
properly represented. 
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A. close scrutiny of the literature reveals that even 
ah initio wave functions fail to represent the experimental 
EPG's correctly. It is well known that, even with carefully 
optimised minimal basis sets, it is not possible to obtain 
reliable values of sensitive one— electron properties such as 
electric field gradients. The goal of most semi-empirical 
formalisms has been to mimic the ab initio SOP MO results. 
Parametrisations and approximations have been done with this 
end in view. It now appears that a different view-point has 
to be taken with regard to the parametrisation procedures in 
the semi— empirical schemes as far as prediction of SPG-'s 
are concerned. Since a large number of accurate experimental 
nuclear quadrupole Coupling constants are becoming available, 
the need for developing a reliable and computationally inex- 
pensive semi-empirical procedure!" with fewer parameters assumes 
significance. 

Since semi-empirical calculations using ’valence only' 
basis sets cannot take into account SternheinBr effects, the 
question of shielding/ antishielding effects in molecules 
with respect to the semi- empirical formalisms needs careful 
scrutiny. It is implicit in all these calculations that 
Sternheimer corrections are either zero or constant. It might 
be possible to develop a theoretical model in ^ich liie atomic- 
like core wave functions are made orthogonal to the semi-empirical 
valence wave functions. Sternheimer shielding calculations can be 
carried out with such core wave functions, as is done in the case 


of atoms. 



APPEMX I 


A Listing of PORTRAM 10 Program for the ETalnation of 
EEG Integrals. is given in this Appendix, fhe Program 
evaluates the Matrix Element of the Diagonal Component of 
the EPG Operator. Atomic Orbitalwise Contributions to EEG 
are also output by this Program. 
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m I'n 30 

TrfMCl.ME.*JC2)GO TO ?9 
TFCMC1.0R.25Ga TO 29 

IFffif?l.Sa.5) JCII, JJ)=4.00»?\t,F1TT3/l5 
TFCNCI.SO*3.0R.MCUE0.4) Qf TT, Jj3»-2.00*ftLFl*T3/t5. 

GO TO 900 

two- IMD THREE CENTRE INTEGRAL 8XALUATIDN 

conttnoe 

CM.rtlLATlON OF INTEGRATinN PARAMETERS AT aUADRAWRR POTNTS 
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Iv the Jn-n K'U^CTIO^*^ RftE EvAj,'M:kT*:n 

* A.JMM with APORTpoiftl'F n RC jv ^\j fc, 
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TF( ARG.GT, 2.05)00 TO 30 
7,=arc;/2.d-) 
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2 -O.2Sl54D-2»T'4nBYX*5 + 0,S3208O-3*TWOB1fX‘‘l 
Sl*St ♦npypf-X) /OSORTCX) 

SUMsSI 

TFC.'I.GT, 1 )m TO 50 , ^ ,v ,:VV' 

BESMOnsSUM 

RFT1JR3 -v 

TF( M3G,G'’’.2.D0)G'3 TO 70 • 

Z?£AR'G/7' 'Di)''' " ' 

SKRIFSs!-:5A«MR-nLOGCZ)+0,5no)*Z '■ ' 

M 0 ■■ ■■ 

COEFFTssSERTES/?. ■■■ 

FACsz .V' f 

MaMtl 


FACBFACMZM*2/0FL0ATfM*(M+n) 
CiFff T wCOSFFTt 0,500 MCI, ODXDr 
TERMeFACMCneFFT 
SERlESaSERieSTTERH 


L0ATCM) + t,»'0/CMfiJ) 



* if • • 


if(ABscre«H/segiES),Gf.i.D*9)sr) to go 
I’rl^^vILTaUD'O/JlRS-iBRlES 

GO TO 80 , : 

APPROXIWATIDM rORMUC#A POR 080E8 1 , X > 2 

XsSilRfl ■ ■ ■■ ' 

S2»f!25ili4p0ff0,234986l9O0»T#O!^XX 

“ - . : 20-l»TWOBXX-2f 0 , l 'j0426RD-ji *T« 1 

■ •O»790l53D-l»TwnayX'‘4+O,32Sgl4r)-2TT.Yn3VK'‘'3 

••9!SR2450-3*T^o6'fx*6 
J2|S2|M|#C,**KJ/0S0RT(X) 

|pfil*l¥?t)G0 ra 90 

feESUdDaWSULTTaRG 

RETURN 

RPCURSIOH PORMOLA POR HIGHPR ORDERS 


Furjf tlaSOM 
KJhJ(?3sRF5nLT 

no t 00 t a? - M 

FllH(t+’l )a2*D0T0Fr..OAT( 1-1 ')1'FilW(T)/&RG ¥ FIMCI"!) 
riiNTTwije ^ : ■ , 

RES 'lODiFUN (« + !)» ARC* *?J 
RETURM 

EiD, ^ ^ . ■ . : 

raNCflOM FACTORCN) : : 

IFCSH,2,3 

RETURN ■ . 

FMol 
GO TO 6 
FMtsI 

no 4 1=1 

FMsF^I^T ' 

COf’ i’l .'lOE 
r ArTOFsK’M 

Pt’TUW''! ' ■ - '■■ 



APPENDIX II 


A List of Conversion Factors for EFGr's and Qnadrupole 
Coupling Constants are given in this Appendix. 



A PPBroiX II 


C omrersios factors 


QGG (in KHz)* = 234. 97x1 0^xQ(in barns )xq (in a.u.) 


Nucleus 

Q( in barns) 

Gonversion factor 
for BBG (in a.u.) 
to QCG (in KHz) 

Conversion factor 
for QGG (in KHz) , 
to EPG (in a.u.) 


2.860x10”^((a) 

672.0 

.1,488x10“^ 

V 

-3.660x10"'^ (b) 

-8599.9 

-1 .1628x10“'^ 


4.099x10“^ (c) 

9631.4 

1.0383x10“^ 

Hr 

8.543x10”^ (c) 

20073.5 

4-982x10“’^ 


1.930x10“^ (d) 

4534.9 

2.2051x10"^ 


-2.630x10“^ (e) 

-6179.7 

-1 .61 8x10”^ 
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